Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1701 (N48E30 L=342)   357    109        0 ( 711/1666- 751/1734) -39.8 E32-E26 N45-N50
      2  691 1619 (N41E30 L=342)   457     67        1 ( 665/1595- 717/1651)  76.0 E31-E26 N39-N43
      3 1072 1481 (N32W02 L= 14)   541    123        1 (1013/1438-1130/1514) -26.6 E02-W06 N29-N34
      4 1341 1417 (N29W21 L= 33)   245     67        1 (1321/1384-1363/1443)  45.0 W19-W23 N26-N30
      5 1594 1322 (N22W40 L= 51)  1915    409        4 (1427/1238-1744/1415) -24.6 W27-W51 N16-N28
      6 1173 1301 (N21W09 L= 21)   913    223        3 (1123/1184-1226/1383) -71.6 W05-W12 N13-N25
      7  637 1343 (N22E27 L=345)   356     64        1 ( 618/1306- 659/1382) -68.2 E29-E25 N19-N24
      8 1337 1261 (N18W20 L= 32)   618    156        0 (1275/1192-1407/1328) -45.0 W15-W24 N14-N22
      9 1004 1036 (N04E01 L= 10)   935    286        1 ( 904/ 968-1134/1136) -28.3 E07-W06 N00-N09
     10  846  919 (S04E11 L=  1)   389    137        1 ( 785/ 903- 924/ 939)   0.0 E14-E05 S04-S02
     11 1154  241 (S54W17 L= 28)   436    111        1 (1105/ 222-1197/ 253) -15.3 W11-W22 S55-S52
  

    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