Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1746 (N71E22 L=105)  1717    139        3 ( 496/1687- 629/1797)   4.6 E36-E10 N65-N76
      2 1230 1629 (N36W35 L=162)   584    219        1 (1130/1598-1327/1668) -48.0 W29-W39 N30-N40
      3 1578 1605 (N23W59 L=186)   728    140        2 (1537/1546-1612/1663)  63.4 W55-W61 N18-N26
      4 1414 1453 (N20W37 L=164)   644    267        2 (1298/1399-1524/1484) -47.5 W30-W42 N13-N25
      5  871 1349 (N29E00 L=127)   389    170        1 ( 816/1308- 951/1398) -60.9 E02-W03 N24-N33
      6  331 1297 (N40E46 L= 81)   582    109        3 ( 295/1253- 356/1333) -50.2 E50-E42 N36-N42
      7  779 1250 (N26E09 L=118)   655    231        0 ( 704/1216- 886/1314) -45.0 E12-E02 N21-N31
      8  712 1024 (N15E18 L=109)   197     99        0 ( 688/ 974- 734/1065)  45.0 E20-E16 N12-N16
      9  508  934 (N14E34 L= 93)   589    141        2 ( 445/ 913- 573/ 971) -26.6 E37-E29 N11-N17
     10  346  906 (N15E47 L= 80)   340     86        1 ( 318/ 874- 370/ 943)  53.1 E50-E44 N13-N17
     11  749  894 (N06E19 L=108)   272     83        1 ( 723/ 847- 774/ 922)  78.7 E20-E17 N03-N08
     12 1291  800 (S12W09 L=137)   482    190        0 (1265/ 716-1331/ 886)  52.1 W06-W13 S17-S08
     13  263  752 (N07E60 L= 67)   795     86        2 ( 230/ 726- 304/ 790) -20.6 E63-E55 N05-N09
     14  878  655 (S10E18 L=109)   235     74        1 ( 842/ 622- 911/ 685)  36.9 E20-E15 S11-S08
     15 1569  613 (S33W25 L=152)   422     59        1 (1540/ 576-1586/ 645)  76.0 W22-W26 S35-S30
     16  710  590 (S10E30 L= 97)   290     91        0 ( 694/ 539- 724/ 629)  71.6 E31-E28 S13-S07
     17  513  513 (S11E47 L= 80)   291     61        1 ( 491/ 486- 533/ 544) -90.0 E48-E45 S12-S08
     18  480  450 (S14E53 L= 74)   868    162        2 ( 416/ 415- 551/ 475) -26.6 E60-E48 S17-S12
     19  515  367 (S21E57 L= 71)   386     70        0 ( 481/ 348- 550/ 389)  11.3 E61-E51 S21-S19
  

    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