Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  746 1829 (N42E51 L= 54)  2321    187        3 ( 700/1760- 808/1907) -16.9 E71-E38 N37-N48
      2  437 1519 (N14E52 L= 53)  1553    229        4 ( 369/1453- 484/1621)  90.0 E56-E48 N08-N20
      3  322 1555 (N14E67 L= 38)   676    109        1 ( 278/1517- 367/1586)  40.6 E70-E63 N10-N16
      4  549 1560 (N19E45 L= 60)   380     84        2 ( 520/1540- 598/1578)  36.9 E47-E41 N17-N20
      5  322 1232 (S05E51 L= 54)  1042    325        0 ( 241/1094- 379/1348)  72.6 E55-E47 S14-N02
      6  740 1076 (S06E18 L= 88)   856    263        0 ( 634/1008- 816/1161)  54.0 E23-E14 S11-N00
      7 1896 1052 (N20W65 L=170)  1488    267        0 (1827/ 934-1948/1145) -20.6 W54-W78 N14-N23
      8 1213 1009 (N01W11 L=116)  1049    244        2 (1167/ 901-1255/1119) -51.3 W06-W15 S04-N06
      9 1721  552 (S13W63 L=168)  2751    295        5 (1595/ 481-1830/ 624)  37.9 W53-W74 S20-S06
     10 1622  364 (S26W70 L=175)   936     91        1 (1581/ 341-1659/ 394)   0.0 W63-W80 S27-S24
  

    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