Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1679 (N27E40 L= 28)   372    117        1 ( 715/1627- 744/1744) -45.0 E45-E37 N23-N31
      2  499 1587 (N16E52 L= 16)   554     95        0 ( 477/1545- 537/1625) -21.8 E56-E46 N13-N17
      3  311 1509 (N07E66 L=  2)   750    112        3 ( 289/1458- 325/1561) -54.5 E71-E63 N03-N10
      4  357 1512 (N08E61 L=  8)   430     46        3 ( 334/1490- 379/1535)  11.3 E63-E58 N06-N09
      5 1457 1434 (N32W16 L= 85)   295     79        2 (1428/1401-1483/1468) -45.0 W13-W18 N30-N34
      6  392  763 (S36E38 L= 31)   697    211        0 ( 323/ 720- 480/ 796)  20.6 E45-E29 S39-S32
      7  664  633 (S37E13 L= 56)  1772    333        3 ( 548/ 550- 785/ 746)  53.6 E21-E06 S45-S26
      8 1721  620 (S07W60 L=128)   572     83        1 (1686/ 587-1754/ 650)  45.0 W56-W62 S09-S04
      9 1682  524 (S13W62 L=131)  1083    236        1 (1589/ 450-1763/ 580)  45.0 W58-W70 S19-S08
     10 1390  439 (S27W44 L=112)   946    249        0 (1300/ 386-1479/ 514)  69.0 W40-W47 S32-S19
     11 1670  468 (S16W66 L=134)   605    127        0 (1625/ 436-1724/ 495)  40.6 W63-W70 S19-S13
     12 1230  369 (S36W38 L=106)   513     76        0 (1204/ 342-1261/ 410) -14.0 W33-W41 S38-S34
  

    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