Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  594 1513 (N34E25 L=268)   809    192        2 ( 512/1445- 664/1564) -35.5 E33-E19 N28-N38
      2  453 1477 (N33E37 L=256)   262     51        0 ( 427/1461- 484/1497) -26.6 E39-E33 N31-N34
      3  747 1106 (N06E16 L=278)  1581    386        3 ( 576/1052- 922/1151)  -5.4 E26-E05 N03-N09
      4 1921 1052 (S08W68 L=  1)  1607    124        5 (1872/ 992-1967/1092) -31.0 W60-W77 S11-S05
      5 1410 1045 (S04W23 L=317)   214     38        0 (1392/1018-1431/1068)  45.0 W22-W24 S04-S02
      6 1843  967 (S12W57 L=351)   933    166        1 (1808/ 900-1882/1026)  45.0 W53-W62 S15-S08
      7  167  958 (N04E62 L=231)   440     46        0 ( 140/ 934- 193/ 977) -31.0 E65-E59 N02-N05
      8 1506  874 (S14W29 L=322)   381    154        0 (1467/ 808-1536/ 942) -63.4 W26-W30 S18-S10
      9 1830  748 (S25W59 L=352)  1094    156        2 (1750/ 714-1882/ 798)  12.1 W49-W65 S27-S22
     10  782  762 (S14E17 L=276)   320    103        1 ( 725/ 745- 839/ 777)  -7.1 E21-E13 S14-S12
     11  885  712 (S18E12 L=282)   363    136        1 ( 861/ 643- 906/ 776)  83.7 E14-E10 S22-S13
     12 1157  681 (S23W05 L=298)   585    145        1 (1105/ 633-1202/ 752)  49.4 W01-W08 S25-S18
     13  418  623 (S18E46 L=247)   485     86        1 ( 398/ 579- 439/ 673)  80.5 E48-E44 S21-S15
     14 1743  493 (S42W60 L=353)   627     35        0 (1713/ 477-1773/ 511)  -5.2 W54-W65 S43-S40
     15  524  413 (S33E47 L=246)  2342    424        2 ( 376/ 313- 681/ 468)   2.4 E81-E33 S38-S29
     16  893  378 (S40E19 L=275)   467    109        0 ( 835/ 345- 946/ 408)   5.7 E24-E14 S42-S38
     17 1153  304 (S50W01 L=294)   329     35        1 (1131/ 281-1179/ 324) -45.0 E01-W03 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