Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1430 1690 (N31W50 L=262)   758    157        0 (1352/1668-1506/1709) -24.4 W44-W55 N27-N34
      2  907 1510 (N35W06 L=217)  1721    636        1 ( 647/1448-1204/1558) -15.8 E14-W25 N27-N39
      3 1492 1350 (N10W36 L=248)   973    236        1 (1431/1236-1537/1439)  85.6 W34-W38 N03-N16
      4 1901 1173 (S10W68 L=279)  1245    131        1 (1868/1131-1961/1219) -21.3 W64-W82 S14-S07
      5  407  874 (N08E42 L=170)  1624    525        2 ( 235/ 708- 567/1007) -43.6 E55-E34 S03-N18
      6  221  824 (N09E60 L=151)   511    129        1 ( 186/ 769- 258/ 876) -74.1 E64-E57 N05-N12
      7 1402  758 (S21W16 L=228)   625    139        0 (1364/ 701-1429/ 827)  45.0 W12-W19 S23-S17
      8 1758  650 (S39W46 L=258)  2758    379        3 (1655/ 527-1852/ 811)  33.7 W35-W56 S45-S30
      9  696  433 (S23E37 L=175)   666    207        1 ( 667/ 339- 732/ 527)  47.5 E42-E31 S29-S17
     10  530  410 (S21E51 L=161)  1031    206        0 ( 500/ 330- 565/ 508)  47.5 E56-E43 S27-S15
     11  825  397 (S29E30 L=181)   344     86        0 ( 800/ 353- 851/ 445)  32.0 E34-E26 S31-S25
     12 1272  302 (S50E04 L=208)  1457    296        2 (1153/ 252-1410/ 362)   0.0 E16-W09 S51-S47
  

    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