Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1566 (N29E35 L=225)   566    151        1 ( 547/1500- 606/1633) -83.7 E36-E32 N24-N33
      2  721 1518 (N25E22 L=238)   445     90        1 ( 676/1484- 773/1544)  18.4 E25-E19 N23-N27
      3 1702 1430 (N25W47 L=307)  1866    300        1 (1591/1301-1785/1537) -49.4 W38-W53 N17-N31
      4  449 1479 (N22E42 L=218)   428    116        0 ( 399/1454- 485/1496)  -7.1 E46-E38 N20-N23
      5  268 1341 (N14E55 L=205)   424     73        1 ( 254/1311- 286/1375) -76.0 E57-E53 N11-N15
      6 1644 1252 (N13W40 L=300)  1340    213        2 (1563/1179-1733/1351) -33.7 W34-W47 N08-N18
      7  752 1301 (N11E18 L=242)   266     45        0 ( 726/1271- 774/1329) -53.1 E19-E16 N09-N13
      8 1234 1262 (N11W12 L=272)   389    129        0 (1165/1238-1304/1288)   6.3 W07-W16 N09-N12
      9 1827 1155 (N09W56 L=316)   380     65        1 (1807/1121-1840/1186)  76.0 W54-W57 N06-N10
     10 1869 1108 (N06W61 L=321)  1225    288        1 (1826/1054-1907/1165) -80.5 W56-W66 N03-N09
     11  158 1103 (S01E65 L=195)   664     74        0 ( 136/1065- 180/1144)  45.0 E67-E61 S02-N01
     12  394 1019 (S06E41 L=219)   974    210        0 ( 355/ 939- 476/1077)  53.1 E44-E34 S10-S02
     13  552  963 (S09E29 L=231)   322     93        0 ( 518/ 930- 594/ 996)  38.7 E31-E26 S11-S07
     14  744  799 (S19E17 L=243)  1116    286        2 ( 651/ 686- 873/ 881)  36.3 E23-E08 S25-S13
     15 1720  778 (S15W49 L=309)   635    154        0 (1660/ 735-1761/ 799)  12.5 W44-W53 S17-S13
     16 1545  601 (S27W39 L=299)  1326    236        1 (1439/ 555-1615/ 697)  48.0 W31-W42 S30-S20
     17  603  543 (S36E29 L=231)   745    214        0 ( 551/ 474- 672/ 603)  39.3 E34-E23 S40-S31
     18 1775  567 (S27W61 L=321)   893    112        1 (1726/ 532-1809/ 592)  33.7 W56-W64 S29-S25
  

    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