Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1529 1671 (N28W59 L=166)  1353    222        1 (1448/1585-1613/1716)   0.0 W46-W73 N26-N30
      2 1651 1509 (N16W57 L=164)   508    104        0 (1604/1483-1694/1535) -11.3 W51-W61 N14-N17
      3 1260 1450 (N25W26 L=133)   679    160        1 (1195/1398-1326/1490)  -5.7 W21-W31 N22-N26
      4 1001 1289 (N23W05 L=112)   686    232        1 ( 960/1190-1065/1366)  48.4 W02-W11 N16-N25
      5  244 1272 (N39E57 L= 50)   756    165        0 ( 231/1205- 263/1359) -68.2 E59-E52 N34-N44
      6 1265 1265 (N14W20 L=127)   295    164        0 (1208/1205-1305/1320) -76.0 W18-W21 N10-N18
      7 1899 1118 (S13W65 L=172)  2865    212        7 (1869/1040-1934/1221)  31.6 W60-W75 S17-S08
      8 1842 1154 (S09W59 L=165)   607     97        2 (1816/1102-1862/1196)  51.3 W54-W61 S11-S06
      9 1799 1024 (S14W51 L=157)   885    157        4 (1755/ 959-1837/1098) -78.7 W48-W54 S18-S08
     10 1768 1001 (S15W47 L=154)   313     50        0 (1753/ 968-1786/1028)  56.3 W45-W48 S16-S13
     11 1695  764 (S26W37 L=144)  2171    455        4 (1599/ 610-1778/ 903)  37.4 W27-W45 S33-S18
     12 1053  856 (S03E02 L=104)   184     75        0 (1018/ 841-1093/ 879) -53.1 E03-W00 S04-N00
     13  912  669 (S10E15 L= 92)   265    136        0 ( 862/ 648- 976/ 705) -45.0 E17-E11 S12-S07
     14 1345  621 (S24W09 L=116)   296    138        0 (1296/ 583-1408/ 654)   7.1 W05-W13 S25-S23
     15 1621  551 (S39W29 L=136)   871    158        0 (1555/ 505-1671/ 613)  15.3 W22-W34 S41-S35
     16 1189  497 (S28E04 L=103)   282     67        0 (1173/ 471-1214/ 533) -90.0 E05-E02 S30-S25
     17  373  447 (S14E64 L= 43)  1836    115        4 ( 343/ 380- 409/ 494)  28.1 E73-E57 S18-S10
     18  586  424 (S18E45 L= 61)   386     79        1 ( 571/ 379- 597/ 467)  56.3 E48-E42 S21-S15
  

    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