Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1073 1803 (N44E20 L=354)   540     69        0 (1032/1782-1114/1822)  15.9 E23-E16 N42-N46
      2 1304 1684 (N41W01 L= 15)   380     43        0 (1276/1661-1325/1703) -18.4 W00-W03 N39-N43
      3  578 1571 (N16E42 L=332)   644    177        0 ( 556/1502- 597/1642) -48.4 E45-E37 N11-N20
      4 1175 1344 (N14W01 L= 15)  1012    299        2 (1122/1221-1261/1480) -52.4 E04-W08 N08-N21
      5 1881 1162 (N24W56 L= 70)   727    121        1 (1844/1125-1908/1225) -26.6 W51-W59 N21-N26
      6  252 1076 (S18E53 L=321)   815    168        0 ( 223/ 995- 269/1154) -60.9 E57-E50 S23-S14
      7  722 1045 (S12E18 L=356)   881    199        0 ( 641/ 964- 782/1143)  72.9 E21-E15 S18-S05
      8 1374  949 (S03W22 L= 36)  1555    407        1 (1259/ 783-1465/1129)  84.8 W18-W25 S14-N08
      9 1020  738 (S23W07 L= 20)   730    153        2 ( 936/ 716-1104/ 768)  15.3 W00-W11 S25-S21
     10 1784  646 (S08W63 L= 77)  1697    209        0 (1705/ 577-1842/ 725)  59.7 W56-W67 S13-S01
  

    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