Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1783 1482 (N12W66 L=291)   869     91        2 (1730/1451-1817/1508) -33.7 W61-W71 N10-N14
      2  724 1364 (N27E12 L=213)   421    126        0 ( 671/1321- 774/1408) -54.5 E14-E09 N23-N31
      3 1652 1299 (N05W45 L=269)   888    172        1 (1575/1252-1725/1341)   8.1 W38-W52 N03-N06
      4 1160 1248 (N12W12 L=237)   406     98        0 (1118/1224-1214/1278) -51.3 W10-W14 N09-N14
      5  427 1242 (N26E35 L=190)   283     93        0 ( 378/1219- 461/1273)   8.1 E38-E31 N24-N27
      6 1358 1181 (N04W22 L=247)   335     60        1 (1323/1149-1388/1210) -53.1 W20-W23 N02-N06
      7  414 1125 (N19E37 L=188)  1404    428        3 ( 256/1067- 598/1172) -20.3 E50-E23 N13-N24
      8  875  878 (S04E11 L=213)   230     79        1 ( 839/ 855- 918/ 898) -45.0 E13-E09 S05-S02
      9  827  748 (S10E17 L=208)   884    306        0 ( 747/ 706- 959/ 823) -17.1 E21-E08 S13-S04
     10  505  670 (S09E39 L=185)   485     89        2 ( 470/ 639- 534/ 716)  26.6 E42-E36 S10-S06
     11 1174  654 (S22W02 L=226)   287     94        1 (1148/ 611-1198/ 700)  53.1 W00-W03 S24-S20
     12 1459  541 (S36W19 L=244)  1641    510        1 (1268/ 405-1665/ 684)  15.9 W02-W37 S41-S30
     13  339  536 (S14E59 L=165)  5559    343       13 ( 253/ 427- 417/ 649)  78.7 E72-E49 S21-S06
     14  269  610 (S08E62 L=162)   447     63        1 ( 244/ 587- 293/ 627)   8.1 E65-E58 S09-S07
     15  647  539 (S19E34 L=190)   293     51        0 ( 624/ 501- 670/ 568)  76.0 E35-E32 S21-S17
     16  867  310 (S39E30 L=194)   821    133        0 ( 809/ 264- 909/ 362)  26.6 E36-E24 S41-S35
     17 1037  261 (S47E20 L=204)   447     71        2 ( 999/ 239-1075/ 281)   0.0 E23-E16 S48-S45
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y