Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1151 1315 (N22W14 L=194)   572    232        1 (1040/1287-1251/1344) -30.3 W07-W19 N18-N25
      2 1757  890 (S17W48 L=228)   316     76        2 (1729/ 871-1780/ 924) -76.0 W44-W50 S18-S14
      3 1810  864 (S20W54 L=234)  1150    167        3 (1744/ 793-1869/ 924)  45.0 W46-W61 S24-S16
      4 1653  876 (S15W38 L=218)   323     83        0 (1612/ 842-1681/ 912)  18.4 W34-W40 S16-S13
      5 1696  839 (S18W42 L=222)   390     98        3 (1673/ 814-1734/ 868) -45.0 W39-W45 S20-S16
      6 1521  583 (S31W26 L=206)  1474    334        1 (1451/ 474-1626/ 720)  42.9 W19-W34 S37-S24
      7 1768  651 (S34W51 L=231)   751    148        1 (1733/ 625-1808/ 685)  -5.7 W46-W56 S36-S29
      8 1022  622 (S18E08 L=172)   491    189        0 ( 952/ 592-1088/ 661) -14.0 E11-E02 S19-S13
      9  510  514 (S17E45 L=135)  1454    364        2 ( 371/ 471- 663/ 546)  -4.6 E58-E32 S18-S13
     10  970  365 (S35E18 L=162)  3101    482       11 ( 831/ 209-1122/ 532)  29.1 E39-E03 S46-S24
  

    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