Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  815 1605 (N36E13 L=184)   407     80        0 ( 770/1588- 857/1627)   0.0 E15-E09 N34-N37
      2 1591 1563 (N30W44 L=241)   448    100        1 (1559/1530-1619/1591)  23.2 W41-W48 N27-N31
      3  745 1302 (N16E16 L=180)   286     99        0 ( 699/1285- 807/1323) -15.9 E19-E12 N14-N17
      4  143 1235 (N15E66 L=130)   996    208        1 ( 111/1202- 178/1271) -26.6 E71-E61 N13-N18
      5 1040 1025 (S02W01 L=197)   598    157        0 (1013/ 950-1060/1118)  90.0 W00-W02 S06-N03
      6  212  958 (S02E58 L=138)   358     58        0 ( 193/ 935- 234/ 987) -36.9 E59-E55 S03-N00
      7 1061  798 (S16W01 L=198)   475    160        0 (1008/ 724-1100/ 870)  60.9 E02-W03 S20-S11
      8  677  516 (S32E27 L=169)   446     77        1 ( 643/ 489- 716/ 550) -36.9 E29-E24 S34-S29
      9  533  488 (S33E40 L=157)  1007    267        1 ( 471/ 428- 617/ 536)  21.8 E46-E31 S36-S29
  

    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