Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  279 1269 (N13E53 L=142)  2383    406        4 ( 193/1150- 338/1442)  77.5 E60-E48 N06-N24
      2 1743 1401 (N21W53 L=248)   587    116        0 (1698/1356-1789/1437) -26.6 W49-W57 N18-N23
      3  686 1315 (N15E21 L=173)   611    178        3 ( 607/1292- 780/1341)   0.0 E26-E15 N13-N16
      4  349 1191 (N08E45 L=149)   585    178        1 ( 304/1112- 391/1251)  53.1 E48-E41 N03-N11
      5  911  841 (S14E07 L=187)   962    294        1 ( 843/ 712- 947/ 977)  69.4 E11-E04 S22-S05
      6 1383  819 (S15W23 L=217)   229     82        0 (1360/ 774-1405/ 849) -59.0 W21-W24 S18-S13
      7 1471  799 (S17W29 L=223)   386    140        1 (1412/ 776-1515/ 827)   8.1 W24-W31 S17-S14
      8  903  471 (S38E09 L=185)   736    216        1 ( 794/ 452-1000/ 503)   7.1 E18-E02 S39-S36
      9 1691  443 (S39W63 L=257)   519     70        0 (1658/ 425-1725/ 462)  -3.8 W56-W71 S40-S37
     10  878  303 (S52E15 L=179)   523    117        0 ( 820/ 286- 936/ 320)  -4.8 E20-E08 S53-S50
  

    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