Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1256 1557 (N28W26 L=211)   296     79        0 (1217/1525-1300/1577) -51.3 W24-W28 N25-N30
      2 1719 1419 (N11W55 L=240)   985    126        2 (1656/1378-1776/1454) -36.9 W50-W59 N07-N14
      3  780 1333 (N24E09 L=175)   232     78        0 ( 743/1320- 819/1348) -11.3 E11-E06 N22-N25
      4 1048 1285 (N16W06 L=191)   865    359        0 ( 908/1240-1208/1331) -17.4 E01-W15 N13-N19
      5  112 1233 (N30E74 L=110)   811     73        2 ( 100/1200- 127/1279) -39.8 E80-E69 N27-N32
      6  327 1131 (N20E44 L=141)  1595    287        5 ( 199/1069- 423/1189) -10.3 E57-E35 N15-N25
      7  119 1040 (N18E69 L=116)  3396    204        9 (  81/ 930- 157/1113)  55.0 E77-E62 N11-N22
      8  180  975 (N13E60 L=124)   581     80        1 ( 156/ 948- 208/1014) -63.4 E63-E57 N11-N15
      9  375  940 (N08E42 L=142)   427     78        1 ( 336/ 920- 405/ 968)   9.5 E45-E39 N06-N09
     10 1156  906 (S08W05 L=190)   473    235        0 (1062/ 821-1246/ 947)  10.3 E01-W10 S10-S05
     11  280  854 (N04E52 L=133)   342     80        0 ( 253/ 818- 301/ 894) -76.0 E53-E50 N02-N06
     12  529  653 (S11E38 L=147)  2231    544        5 ( 428/ 478- 587/ 891)  47.7 E52-E29 S20-N02
     13  844  716 (S13E16 L=168)   978    317        2 ( 752/ 635- 941/ 812)  19.7 E22-E08 S17-S08
     14 1220  729 (S19W06 L=191)   345    124        1 (1196/ 656-1244/ 767)  66.8 W04-W08 S23-S16
     15 1738  586 (S40W48 L=232)  1889    365        2 (1599/ 482-1876/ 690)  15.7 W33-W65 S44-S35
     16  300  609 (S10E59 L=126)   661    105        2 ( 270/ 572- 329/ 659)  20.0 E64-E53 S11-S07
     17  624  587 (S17E34 L=151)   498     91        2 ( 599/ 534- 649/ 642)  45.0 E36-E30 S19-S13
     18  680  456 (S26E35 L=150)  1097    291        1 ( 549/ 424- 815/ 497)  -2.7 E45-E24 S27-S23
     19 1395  367 (S48W14 L=199)   554     65        1 (1365/ 337-1427/ 394)   9.5 W11-W17 S50-S46
  

    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