Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  475 1537 (N29E49 L=213)  1892    330        4 ( 434/1423- 559/1649) -35.5 E61-E37 N21-N36
      2 1449 1520 (N38W29 L=290)   705    155        0 (1392/1461-1502/1567) -41.2 W24-W32 N33-N41
      3 1157 1202 (N14W07 L=269)   240    121        0 (1100/1168-1198/1219) -26.6 W03-W09 N11-N14
      4 1255  951 (S01W15 L=277)   598    129        1 (1189/ 917-1309/ 983)  23.2 W10-W18 S03-N01
      5 1687  938 (N02W46 L=308)   543     90        1 (1655/ 909-1723/ 969)  26.6 W43-W49 N00-N03
      6  276  853 (S16E53 L=208)   731     89        2 ( 236/ 816- 308/ 888) -26.6 E57-E49 S18-S13
      7  260  770 (S22E57 L=205)  1106    162        6 ( 238/ 687- 287/ 823)  56.3 E62-E53 S27-S18
      8 1225  738 (S14W15 L=277)   978    320        4 (1094/ 672-1352/ 810) -15.5 W06-W24 S17-S10
      9 1006  726 (S17W01 L=263)   327    122        0 ( 954/ 701-1057/ 740)  21.8 E01-W04 S18-S15
     10  461  675 (S26E37 L=224)   622    150        0 ( 399/ 658- 532/ 696)  18.4 E43-E31 S27-S23
     11 1799  610 (S19W69 L=331)   609     49        0 (1781/ 584-1815/ 639) -76.0 W66-W71 S20-S16
     12 1107  471 (S33W12 L=274)   244     60        1 (1068/ 449-1132/ 489) -11.3 W08-W13 S35-S32
  

    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