Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1350 1498 (N25W35 L=269)  1175    330        2 (1211/1449-1499/1541) -15.3 W23-W46 N21-N27
      2  802 1484 (N40E02 L=232)   685    149        2 ( 734/1457- 876/1511) -11.3 E07-W03 N37-N42
      3  595 1466 (N44E20 L=214)   243     50        3 ( 575/1439- 613/1490) -63.4 E21-E18 N42-N46
      4 1749 1463 (N11W65 L=299)   422     79        2 (1704/1449-1781/1481) -11.3 W59-W69 N09-N12
      5 1480 1414 (N17W39 L=273)   312     63        0 (1458/1389-1504/1442) -76.0 W38-W41 N14-N18
      6  900 1220 (N21E03 L=231)   805    368        1 ( 751/1130-1036/1281) -45.0 E11-W03 N12-N27
      7  378 1107 (N26E42 L=192)  1764    392        4 ( 226/1009- 510/1235) -35.1 E59-E32 N17-N36
      8  383  863 (N11E45 L=189)   318    140        0 ( 337/ 833- 447/ 878)   0.0 E49-E39 N09-N12
      9  720  752 (S02E24 L=210)   440    138        0 ( 679/ 703- 745/ 818)  29.7 E28-E21 S03-N01
     10 1560  748 (S23W26 L=260)   519    199        1 (1484/ 682-1621/ 808)  14.0 W19-W31 S25-S21
     11 1404  711 (S21W15 L=249)   757    190        1 (1369/ 619-1430/ 790)  52.1 W10-W18 S25-S15
     12 1744  534 (S47W47 L=281)   936     81        1 (1700/ 507-1782/ 555) -23.6 W39-W55 S49-S42
     13 1577  413 (S50W26 L=260)   699     70        2 (1546/ 378-1605/ 448)  71.6 W22-W30 S52-S46
     14  447  374 (S20E63 L=171)   605     71        1 ( 423/ 350- 477/ 404)  14.0 E69-E57 S22-S18
     15 1495  370 (S50W17 L=251)   552     60        2 (1473/ 341-1514/ 404)  59.0 W15-W19 S52-S46
     16 1221  333 (S42E07 L=227)   537     65        1 (1190/ 302-1250/ 358)  38.7 E09-E04 S44-S39
     17  739  314 (S30E42 L=192)   332     46        1 ( 714/ 296- 764/ 331) -56.3 E43-E39 S31-S28
     18  803  255 (S36E42 L=192)   694    106        2 ( 768/ 215- 826/ 294)  20.0 E48-E37 S39-S33
  

    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