Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1381 1785 (N54W52 L= 43)  1474    217        3 (1300/1753-1453/1815)   3.6 W38-W70 N51-N56
      2  450 1532 (N41E47 L=304)   781    120        1 ( 397/1504- 507/1553)   0.0 E54-E41 N38-N42
      3  814 1450 (N33E12 L=339)   618    222        0 ( 783/1349- 846/1532) -81.3 E14-E09 N26-N39
      4  557 1496 (N38E34 L=317)   318     49        2 ( 533/1465- 584/1519) -63.4 E37-E32 N35-N39
      5 1048 1448 (N32W05 L=356)   999    211        3 ( 969/1413-1151/1498) -20.0 W00-W11 N28-N35
      6 1070 1236 (N18W04 L=356)   689    220        2 (1010/1154-1152/1300)  38.7 W00-W10 N12-N21
      7 1618 1153 (N08W41 L= 32)   512     90        3 (1581/1109-1646/1184) -90.0 W38-W42 N05-N09
      8 1525 1092 (N05W33 L= 24)   604    133        0 (1448/1075-1586/1111) -11.3 W27-W37 N03-N06
      9 1770  788 (S17W53 L= 45)   437     54        0 (1750/ 768-1796/ 818) -38.7 W50-W56 S18-S14
     10 1223  612 (S23W10 L=  2)   227    104        0 (1185/ 578-1263/ 650)  33.7 W07-W13 S24-S20
     11  493  499 (S27E44 L=307)  1347    362        0 ( 426/ 408- 590/ 563) -35.0 E57-E36 S33-S21
  

    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