Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  830 1855 (N41E50 L=  8)  1913    185        4 ( 761/1805- 905/1924) -10.2 E79-E39 N37-N46
      2  511 1652 (N18E56 L=  2)  1511    258        1 ( 444/1545- 589/1749) -76.0 E60-E52 N09-N25
      3  805 1680 (N27E36 L= 22)   590    109        0 ( 766/1639- 860/1714)  -5.7 E40-E30 N25-N28
      4 1575 1651 (N54W22 L= 80)   406     51        0 (1550/1627-1593/1669) -11.3 W19-W24 N51-N56
      5 1401 1536 (N36W08 L= 66)   452    114        0 (1360/1484-1442/1587) -26.6 W03-W12 N33-N38
      6 1200 1536 (N29E05 L= 53)   287     86        0 (1166/1493-1218/1577) -81.9 E06-E03 N25-N32
      7  314 1519 (N05E66 L=351)   682     78        0 ( 285/1493- 351/1549)  -9.5 E72-E60 N04-N06
      8  622 1500 (N11E38 L= 20)   250     64        0 ( 594/1469- 649/1524)  76.0 E39-E37 N08-N12
      9 1013 1475 (N20E14 L= 44)   314     94        1 ( 970/1457-1061/1492)   9.5 E16-E10 N18-N20
     10 1638 1304 (N29W31 L= 89)   267     75        0 (1606/1267-1668/1326)  -9.5 W28-W34 N27-N30
     11  581 1282 (S03E33 L= 25)   409    126        1 ( 517/1233- 627/1321)  56.3 E35-E30 S06-N00
     12 1567 1205 (N21W27 L= 85)   236     45        1 (1548/1182-1580/1235) -56.3 W26-W28 N18-N22
     13  165 1034 (S26E66 L=352)   896    117        3 ( 121/ 991- 206/1081)  30.3 E74-E60 S29-S22
     14  285 1036 (S24E51 L=  7)   385     50        1 ( 262/1013- 308/1061) -71.6 E53-E48 S24-S21
     15  441 1001 (S22E36 L= 22)   321     76        1 ( 418/ 963- 461/1037) -76.0 E37-E33 S23-S19
     16 1896  947 (N18W65 L=123)  1295    160        1 (1853/ 904-1946/1017) -13.0 W57-W77 N14-N20
     17 1177  971 (S04W10 L= 68)   217     68        1 (1155/ 939-1198/1008) -76.0 W09-W10 S05-S01
     18 1715  655 (S04W57 L=115)   567    112        1 (1674/ 620-1750/ 677)  45.0 W53-W60 S07-S02
     19  894  585 (S34W06 L= 64)   871    201        2 ( 840/ 539- 970/ 631)   5.2 W00-W11 S38-S29
  

    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