Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1180 1524 (N34W02 L= 12)  1572    344        1 (1040/1426-1341/1613) -16.9 E09-W14 N28-N37
      2  779 1535 (N28E26 L=345)   510     91        2 ( 756/1496- 818/1579) -71.6 E27-E22 N24-N30
      3 1621 1470 (N38W38 L= 49)   336    108        0 (1579/1453-1665/1484)   0.0 W34-W42 N36-N38
      4 1400 1403 (N29W19 L= 30)   230     29        1 (1384/1379-1418/1426) -71.6 W18-W20 N27-N30
      5  764 1332 (N14E21 L=349)   290    150        0 ( 711/1286- 818/1393) -29.1 E25-E16 N12-N17
      6  309 1313 (N07E55 L=315)   443     67        0 ( 286/1282- 340/1340)  36.9 E56-E52 N04-N08
      7 1733 1212 (N22W47 L= 57)  1250    268        1 (1651/1138-1814/1276) -20.6 W39-W55 N18-N24
      8 1149 1185 (N11W05 L= 15)   687    175        3 (1068/1098-1198/1266)  71.6 W01-W07 N04-N16
  

    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