Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  504 1564 (N15E47 L=224)  1150    247        4 ( 463/1453- 538/1665) -68.2 E54-E42 N07-N22
      2  454 1592 (N17E53 L=218)   547     62        4 ( 427/1566- 483/1620)  -8.1 E56-E49 N14-N18
      3 1314 1541 (N31W07 L=277)   489    151        1 (1276/1479-1357/1608)  79.7 W05-W08 N25-N36
      4 1484 1465 (N30W20 L=291)  1291    375        0 (1381/1317-1577/1598) -39.5 W10-W28 N22-N37
      5  585 1533 (N15E39 L=231)   545     76        2 ( 549/1495- 619/1570) -78.7 E42-E37 N12-N17
      6  896 1459 (N16E17 L=254)   696    176        1 ( 836/1428- 981/1508)  -5.7 E21-E11 N14-N17
      7 1552 1104 (N09W30 L=300)  1475    318        5 (1518/ 955-1600/1241) -76.8 W27-W33 N01-N18
      8 1679 1107 (N13W39 L=310)   386    115        0 (1659/1049-1693/1157) -59.0 W37-W41 N09-N16
      9 1811  959 (N08W54 L=325)  1194    190        0 (1734/ 907-1861/1021)  39.3 W47-W58 N04-N13
     10 1606  923 (N01W37 L=308)   304     78        1 (1570/ 898-1644/ 944)   0.0 W34-W40 N00-N01
     11 1360  744 (S15W27 L=297)   503    246        0 (1264/ 699-1474/ 793)  36.0 W21-W33 S19-S10
     12 1565  605 (S18W46 L=316)   699    115        1 (1514/ 568-1614/ 646)  63.4 W42-W47 S20-S14
     13 1710  590 (S14W59 L=330)   496     60        1 (1681/ 558-1737/ 620) -56.3 W55-W61 S15-S12
     14 1707  521 (S18W64 L=334)  1102    175        3 (1630/ 484-1780/ 560)  36.0 W57-W71 S21-S13
     15 1365  466 (S31W39 L=309)   711    141        0 (1307/ 433-1426/ 515)  63.4 W35-W41 S34-S26
     16 1232  426 (S37W32 L=302)   258     12        1 (1213/ 404-1250/ 448) -45.0 W29-W33 S38-S35
     17 1457  304 (S38W60 L=330)   520     45        0 (1430/ 286-1487/ 327)  -5.7 W54-W64 S38-S36
  

    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