Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1024 1725 (N35E22 L=343)  1257    201        2 ( 961/1643-1071/1798) -29.4 E31-E15 N30-N39
      2  653 1645 (N20E42 L=323)   937    220        1 ( 570/1597- 748/1675)  31.0 E45-E35 N15-N23
      3  469 1534 (N09E50 L=316)   769    155        0 ( 445/1466- 495/1602) -48.8 E53-E45 N05-N13
      4 1421 1481 (N31W12 L= 17)  1200    241        0 (1336/1398-1493/1537) -17.1 W05-W18 N26-N33
      5  333 1453 (N02E58 L=307)   382     93        0 ( 317/1410- 347/1494) -33.7 E61-E55 N00-N04
      6 1067 1351 (N12E06 L=359)   594    156        0 (1026/1289-1131/1408)  71.6 E07-E03 N07-N16
      7 1287 1249 (N12W09 L= 14)   734    261        1 (1207/1149-1361/1350) -26.6 W01-W15 N08-N16
      8  219 1240 (S12E62 L=303)   405     91        0 ( 182/1225- 259/1258)  18.4 E66-E57 S13-S10
      9  248 1010 (S25E54 L=311)   332    100        0 ( 213/ 979- 270/1050)  39.8 E58-E52 S27-S22
     10  341  831 (S35E43 L=322)   808    200        0 ( 258/ 780- 407/ 894)  32.0 E53-E37 S39-S29
     11 1534  670 (S12W41 L= 47)   600    174        0 (1457/ 615-1605/ 708)  56.3 W38-W44 S17-S07
     12 1068  631 (S28W14 L= 20)   858    197        2 ( 989/ 595-1145/ 703)  61.4 W10-W17 S32-S21
     13 1683  589 (S11W57 L= 62)   460     65        1 (1664/ 556-1711/ 616) -76.0 W54-W58 S13-S09
     14  349  566 (S53E47 L=319)   347     29        0 ( 328/ 547- 369/ 587)  20.6 E50-E42 S54-S50
     15 1621  436 (S21W63 L= 68)  2675    130        7 (1566/ 362-1703/ 487)  -8.5 W54-W74 S24-S16
     16 1374  444 (S29W42 L= 48)   810    160        0 (1290/ 426-1453/ 462)  22.6 W36-W48 S31-S26
     17 1161  348 (S42W35 L= 41)   362     61        1 (1137/ 320-1180/ 383) -38.7 W33-W38 S43-S39
  

    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