Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1691 1488 (N15W60 L=273)   999    170        2 (1651/1401-1726/1559)  63.4 W54-W64 N09-N19
      2  324 1421 (N46E50 L=163)   390     85        2 ( 292/1402- 365/1437) -10.3 E55-E44 N44-N47
      3 1458 1265 (N10W32 L=245)   223     88        0 (1439/1231-1483/1308) -90.0 W30-W32 N07-N12
      4 1414 1037 (S01W23 L=236)  1083    388        2 (1307/ 863-1518/1138)  25.3 W12-W31 S08-N03
      5 1786  995 (S14W49 L=263)   475     46        1 (1759/ 967-1812/1021)  21.8 W46-W52 S15-S12
      6  909  755 (S06E13 L=201)   343    106        2 ( 897/ 713- 929/ 817)  81.9 E14-E11 S08-S01
      7  970  755 (S07E09 L=204)   199    109        0 ( 937/ 730-1012/ 799) -71.6 E10-E07 S09-S03
      8 1158  766 (S11W02 L=215)   317    100        1 (1107/ 744-1206/ 786)   0.0 E01-W04 S11-S09
      9 1124  703 (S13E02 L=212)   231    112        1 (1093/ 676-1159/ 729)  18.4 E03-W00 S14-S12
     10  543  613 (S07E39 L=174)   383     78        1 ( 511/ 578- 576/ 648)  78.7 E41-E37 S09-S04
     11  620  564 (S11E36 L=178)   651    167        0 ( 589/ 505- 668/ 630) -83.7 E37-E34 S15-S06
     12  529  520 (S13E44 L=169)  1138    196        6 ( 473/ 479- 605/ 585) -49.4 E49-E38 S15-S07
     13  952  528 (S20E16 L=197)   258     95        0 ( 925/ 493- 983/ 577)  45.0 E18-E13 S22-S17
     14  512  437 (S18E50 L=164)   961    160        2 ( 442/ 393- 587/ 476) -36.9 E56-E45 S21-S14
     15  737  337 (S29E39 L=175)  1080    288        2 ( 653/ 253- 845/ 442)  21.0 E51-E25 S34-S23
     16  507  367 (S22E56 L=158)   875    143        3 ( 479/ 313- 542/ 417)  63.4 E59-E50 S26-S18
     17  463  337 (S24E64 L=150)  1933    220        5 ( 396/ 234- 523/ 413)  41.6 E76-E57 S34-S18
  

    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