Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  668 1809 (N46E44 L=262)  1456    154        2 ( 613/1776- 744/1843) -12.8 E56-E34 N42-N49
      2 1095 1703 (N40E02 L=304)  1237    418        0 ( 896/1677-1262/1738)   6.3 E17-W10 N36-N41
      3  625 1628 (N30E36 L=271)   609    148        1 ( 584/1567- 649/1677) -45.0 E41-E31 N25-N33
      4 1522 1648 (N42W34 L=340)   454     89        0 (1476/1628-1564/1668)  -6.3 W29-W38 N40-N43
      5  781 1598 (N29E23 L=283)   518    100        3 ( 762/1548- 807/1655) -63.4 E25-E20 N25-N33
      6 1764 1546 (N39W61 L=  7)  1177    143        2 (1704/1489-1808/1586) -24.4 W53-W67 N34-N40
      7  307 1469 (N18E58 L=248)  1562    159        3 ( 270/1406- 362/1538) -27.9 E67-E50 N13-N22
      8 1099 1352 (N15W02 L=308)   531    228        0 (1055/1260-1148/1464) -63.4 E01-W05 N09-N22
      9 1557 1310 (N18W32 L=338)   359    148        0 (1497/1285-1629/1345) -11.3 W27-W37 N17-N19
     10 1179 1221 (N08W07 L=314)   188     62        0 (1146/1207-1211/1233)  14.0 W05-W09 N07-N08
     11 1483 1181 (N09W27 L=333)   326    115        0 (1450/1127-1518/1227) -51.3 W24-W29 N06-N11
     12 1829 1142 (N11W56 L=  2)   420     76        0 (1814/1111-1847/1182)  76.0 W54-W58 N09-N13
     13 1486 1047 (N01W28 L=334)   537    115        0 (1442/1008-1542/1074)  26.6 W25-W31 S01-N03
     14 1789  943 (S01W53 L=360)  1104    161        1 (1727/ 903-1864/ 985)  21.0 W48-W61 S04-N02
     15  441  626 (S34E41 L=265)   682    174        0 ( 420/ 536- 468/ 703)  59.7 E45-E36 S40-S28
     16 1038  506 (S37W07 L=313)  2855    481        5 ( 864/ 404-1208/ 664)  47.9 E05-W15 S46-S25
     17 1652  641 (S21W49 L=355)   398     97        1 (1606/ 617-1703/ 662)  26.6 W45-W53 S22-S18
     18  650  595 (S35E23 L=283)   568    142        0 ( 588/ 574- 722/ 613)  10.3 E28-E17 S36-S33
     19 1056  373 (S47W11 L=317)   250     60        0 (1029/ 357-1091/ 388)  26.6 W08-W13 S48-S45
  

    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