Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  564 1738 (N24E55 L=303)   620    186        0 ( 519/1711- 609/1765)  -4.8 E61-E49 N22-N25
      2 1155 1206 (N07W02 L=  0)   371    115        1 (1116/1158-1177/1265)  90.0 W01-W03 N03-N11
      3  276 1183 (S15E53 L=305)   544     88        2 ( 260/1135- 292/1234) -71.6 E55-E50 S18-S12
      4  741 1155 (S06E19 L=339)   368     76        1 ( 709/1113- 780/1185)  53.1 E20-E17 S09-S04
      5 1787 1092 (N21W47 L= 45)  1127    233        0 (1729/ 979-1837/1174) -26.6 W39-W53 N16-N23
      6  265 1058 (S23E51 L=307)   911    230        0 ( 201/1027- 343/1096)   0.0 E58-E43 S24-S21
      7  900  760 (S24W00 L=359)   362     73        0 ( 850/ 742- 937/ 784)  51.3 E02-W02 S26-S21
      8  985  588 (S32W11 L= 10)  1313    633        4 ( 802/ 503-1177/ 688)   2.0 E03-W26 S33-S28
      9 1659  436 (S17W64 L= 63)  1390    218        0 (1580/ 361-1728/ 527)  77.0 W61-W66 S23-S10
     10 1316  323 (S36W46 L= 44)   378     77        1 (1280/ 309-1364/ 338)  45.0 W44-W48 S37-S33
  

    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