Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  379 1580 (N10E62 L= 17)  2617    485        9 ( 313/1467- 446/1664)  85.9 E68-E54 N02-N16
      2  808 1615 (N21E31 L= 48)   312    116        0 ( 766/1573- 840/1654) -20.6 E35-E27 N19-N23
      3 1245 1369 (N19W03 L= 81)   857    284        0 (1218/1243-1275/1500) -68.2 W00-W07 N11-N27
      4 1592 1364 (N30W25 L=104)   296     50        1 (1568/1342-1617/1390)   0.0 W23-W27 N29-N31
      5  242 1316 (S08E62 L= 17)   923    201        1 ( 212/1284- 288/1349)  -5.2 E66-E55 S10-S06
      6  741 1278 (N00E23 L= 56)   536    170        0 ( 659/1238- 830/1300)   5.2 E27-E16 S01-N01
      7 1439 1213 (N16W18 L= 97)   474    131        0 (1418/1149-1460/1279) -56.3 W15-W20 N12-N18
      8  170  958 (S31E59 L= 19)  1608    416        2 (  90/ 934- 252/ 979)  13.2 E67-E50 S32-S28
      9 1255  923 (S06W15 L= 94)   469     95        1 (1232/ 874-1275/ 978) -80.5 W13-W16 S08-S02
     10  672  874 (S24E17 L= 62)   735    151        1 ( 612/ 806- 727/ 940)  70.0 E19-E14 S29-S18
     11 1831  857 (N10W58 L=137)   501     92        1 (1816/ 816-1848/ 897) -53.1 W55-W61 N07-N11
     12  429  757 (S38E33 L= 46)  1367    295        3 ( 324/ 722- 546/ 807)  19.3 E42-E22 S42-S33
     13 1827  685 (N01W67 L=145)  1171    267        0 (1787/ 658-1884/ 743)  15.9 W60-W76 S01-N03
     14  918  605 (S33W07 L= 85)   931    258        2 ( 841/ 532-1019/ 716)  79.4 W02-W09 S39-S23
     15 1261  599 (S23W27 L=106)   504    132        0 (1190/ 565-1319/ 648)  58.0 W23-W28 S26-S18
     16 1551  587 (S14W46 L=125)   348    107        1 (1507/ 572-1597/ 603)  23.2 W42-W49 S15-S12
     17 1280  463 (S30W35 L=114)   342     80        2 (1241/ 446-1327/ 479)  31.0 W32-W37 S30-S27
     18 1426  313 (S32W54 L=133)   385     66        1 (1397/ 300-1457/ 333)  26.6 W52-W56 S34-S31
  

    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