Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  863 1645 (N28E24 L=100)  1514    352        7 ( 796/1492- 957/1794) -42.3 E37-E14 N18-N38
      2  515 1741 (N30E58 L= 66)  1176    149        1 ( 475/1690- 566/1789) -16.3 E72-E48 N26-N33
      3 1244 1612 (N34W02 L=126)   261     64        0 (1221/1592-1272/1639) -33.7 W00-W03 N32-N35
      4  732 1371 (N08E24 L=101)  1050    381        3 ( 704/1218- 756/1577) -70.0 E29-E20 S01-N21
      5 1502 1415 (N27W23 L=147)   406     95        0 (1485/1372-1522/1464)  90.0 W21-W24 N23-N30
      6 1074 1431 (N18E05 L=120)   222     60        1 (1042/1417-1106/1454)   0.0 E07-E03 N16-N19
      7 1107 1394 (N16E02 L=122)   547    153        2 (1028/1369-1177/1416)   6.3 E07-W02 N14-N18
      8 1652 1373 (N28W35 L=160)   437    109        0 (1628/1314-1676/1414) -56.3 W33-W37 N24-N30
      9 1050  907 (S13W04 L=128)   591    182        1 ( 990/ 884-1148/ 938)  24.0 W00-W09 S15-S10
     10  624  823 (S26E22 L=103)   328     95        0 ( 572/ 789- 660/ 849)  -8.1 E25-E18 S27-S24
     11  573  689 (S35E24 L=100)   326     74        1 ( 542/ 669- 610/ 719) -11.3 E26-E21 S36-S33
     12 1422  498 (S28W40 L=165)   734    187        3 (1344/ 442-1475/ 562) -14.9 W32-W47 S30-S26
     13 1247  522 (S31W27 L=151)   336     85        0 (1203/ 509-1291/ 535)  18.4 W23-W29 S32-S29
     14  740  487 (S46E09 L=116)   803    150        2 ( 682/ 464- 831/ 517)  14.0 E13-E01 S47-S42
     15 1421  368 (S36W49 L=173)   519     79        1 (1389/ 342-1463/ 398)  68.2 W45-W50 S38-S33
     16  746  233 (S68E03 L=122)   858    130        1 ( 684/ 203- 810/ 258)  -2.7 E13-W08 S70-S65
  

    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