Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1096 1770 (N52E10 L=162)   643    127        0 (1034/1745-1158/1790)  31.0 E14-E04 N48-N54
      2  540 1755 (N39E62 L=110)   530     87        0 ( 500/1736- 578/1769)  23.2 E66-E56 N37-N40
      3  752 1561 (N29E28 L=144)  2102    407        2 ( 590/1460- 919/1640)  -9.5 E43-E13 N25-N32
      4 1624 1560 (N44W43 L=215)   798     78        0 (1572/1535-1667/1583)  17.1 W36-W49 N41-N47
      5  933 1292 (N15E10 L=162)   403    144        0 ( 882/1224- 972/1338) -35.5 E13-E06 N11-N16
      6  308 1295 (N06E54 L=118)   546     77        1 ( 289/1255- 337/1336) -51.3 E56-E50 N03-N08
      7 1735 1169 (N19W47 L=219)   289     46        0 (1711/1145-1759/1190) -26.6 W45-W49 N17-N19
      8 1896 1020 (N12W67 L=239)  5063    264        3 (1861/ 898-1940/1122) -77.9 W62-W76 N04-N18
      9 1803  936 (N06W56 L=228)   788    125        0 (1788/ 865-1819/1000) -63.4 W53-W58 N01-N09
     10  948  842 (S12E02 L=170)   370    132        0 ( 892/ 811- 999/ 888)  50.2 E05-W00 S14-S08
     11 1030  666 (S22W06 L=178)  1781    530        2 ( 793/ 620-1186/ 790)  36.0 E09-W14 S27-S11
     12  984  292 (S50W12 L=184)   381     92        0 ( 943/ 281-1033/ 311)   0.0 W08-W17 S51-S48
     13 1148  216 (S54W33 L=205)   524     55        0 (1121/ 196-1167/ 246) -59.0 W29-W37 S55-S50
  

    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