Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  591 1500 (N33E25 L=269)   668    151        0 ( 520/1443- 649/1544) -36.0 E31-E20 N28-N37
      2  432 1471 (N33E39 L=256)   334     40        0 ( 407/1455- 467/1492) -26.6 E41-E35 N31-N34
      3  722 1099 (N06E17 L=278)  1702    397        4 ( 536/1050- 900/1144)  -2.5 E29-E06 N02-N09
      4 1892 1027 (S09W63 L=358)  3312    374        6 (1796/ 897-1962/1137)  37.4 W52-W75 S15-S02
      5  152  952 (N04E64 L=230)   436     50        1 ( 133/ 925- 180/ 974) -26.6 E67-E60 N02-N05
      6 1484  871 (S14W27 L=322)   488    159        0 (1444/ 802-1520/ 943) -66.0 W24-W29 S18-S09
      7 1386  817 (S17W20 L=315)   196     75        0 (1349/ 798-1419/ 833)  11.3 W17-W22 S17-S15
      8 1791  737 (S26W54 L=349)   856    163        0 (1713/ 712-1855/ 785)   7.6 W45-W61 S27-S23
      9  764  755 (S14E19 L=276)   370    114        1 ( 720/ 725- 820/ 779)  -9.5 E21-E15 S15-S13
     10  862  708 (S18E13 L=281)   441    128        1 ( 838/ 639- 887/ 771)  69.4 E15-E11 S22-S14
     11 1125  663 (S24W03 L=297)   546    117        1 (1080/ 626-1168/ 714)  45.0 W00-W05 S25-S20
     12  401  614 (S19E48 L=247)   628    112        1 ( 378/ 569- 425/ 671)  50.2 E51-E45 S21-S15
     13 1735  491 (S42W58 L=353)   777     57        1 (1704/ 472-1770/ 509)  -4.8 W53-W65 S43-S40
     14  478  428 (S31E50 L=244)   731    120        2 ( 438/ 377- 525/ 477)  28.6 E56-E44 S34-S28
     15  609  440 (S32E38 L=256)   659    175        0 ( 527/ 388- 666/ 466) -33.7 E44-E34 S36-S29
     16  950  399 (S39E14 L=281)   957    279        0 ( 826/ 341-1074/ 465)  12.8 E25-E03 S42-S36
     17 1137  302 (S50E01 L=294)   316     44        1 (1110/ 281-1162/ 322) -26.6 E03-W01 S51-S47
  

    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