Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  657 1709 (N29E46 L=137)  1564    306        3 ( 542/1659- 745/1742)   6.0 E57-E38 N26-N30
      2  409 1617 (N17E61 L=121)  2045    193        3 ( 339/1559- 469/1678)  52.1 E67-E56 N13-N22
      3 1147 1572 (N32E08 L=175)   307     93        0 (1102/1530-1184/1598) -15.9 E10-E03 N30-N33
      4 1670  739 (S03W49 L=231)   741    112        3 (1639/ 688-1695/ 788) -51.3 W46-W51 S05-N00
      5 1345  708 (S13W27 L=210)   439    127        0 (1299/ 651-1404/ 765)  77.5 W25-W29 S17-S08
      6  732  544 (S40E06 L=176)   838    199        0 ( 641/ 526- 833/ 563)  15.9 E13-W01 S42-S37
      7  329  518 (S54E48 L=135)   484    112        0 ( 293/ 502- 362/ 534)   6.3 E52-E43 S55-S53
      8 1255  398 (S34W35 L=218)   378    104        1 (1212/ 373-1310/ 428)   8.1 W31-W38 S35-S31
      9 1589  382 (S25W61 L=243)  1547    411        1 (1499/ 338-1656/ 413)   6.0 W52-W71 S27-S22
  

    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