Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1621 1251 (N15W42 L=336)   383    128        0 (1565/1215-1670/1279) -35.5 W38-W45 N12-N17
      2 1324 1146 (N11W19 L=313)   312     58        1 (1304/1116-1350/1172)  63.4 W17-W21 N08-N12
      3 1417 1099 (N07W25 L=319)   251     92        0 (1387/1074-1459/1123)  31.0 W23-W28 N05-N08
      4  401 1047 (N06E42 L=252)   345    100        0 ( 365/1007- 432/1083) -39.8 E45-E39 N03-N08
      5 1868  858 (S11W66 L=360)  1291    177        2 (1821/ 781-1902/ 921) -37.6 W58-W74 S16-S06
      6  455  850 (S06E38 L=255)   912    223        0 ( 356/ 787- 532/ 885) -21.0 E46-E32 S09-S03
      7  258  659 (S19E62 L=231)   639     96        2 ( 242/ 617- 279/ 723)  66.8 E65-E59 S21-S14
      8  531  605 (S22E36 L=257)   296     93        0 ( 512/ 561- 549/ 647)  71.6 E38-E34 S24-S18
      9  763  519 (S28E20 L=274)   466    114        1 ( 721/ 473- 799/ 563)  41.2 E24-E16 S31-S24
     10  498  516 (S28E42 L=252)   312     91        0 ( 465/ 499- 533/ 548)  15.9 E44-E37 S29-S25
     11 1632  501 (S33W49 L=342)   717     98        0 (1596/ 473-1677/ 540) -14.0 W43-W55 S36-S30
  

    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