Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1583 1575 (N23W56 L=106)   707    107        0 (1537/1532-1625/1611) -71.6 W53-W57 N20-N26
      2 1243 1546 (N32W29 L= 80)   521     78        2 (1210/1506-1289/1579) -59.0 W26-W31 N28-N33
      3 1637 1296 (N07W45 L= 96)  1916    514        3 (1587/1076-1695/1485)  61.4 W39-W53 S05-N17
      4 1764 1380 (N07W60 L=111)  1541    159        4 (1698/1338-1843/1427) -31.6 W54-W68 N02-N11
      5  522 1348 (N37E28 L= 23)   304     71        1 ( 488/1318- 548/1372) -45.0 E30-E25 N34-N39
      6  769 1266 (N27E10 L= 40)   307    105        3 ( 741/1227- 807/1298) -63.4 E11-E08 N23-N29
      7 1114 1197 (N15W09 L= 60)   519    265        1 (1054/1105-1183/1260)  31.0 W04-W14 N09-N17
      8  276 1138 (N28E52 L=359)   305     67        0 ( 254/1102- 294/1163) -38.7 E54-E49 N25-N29
      9  467 1071 (N21E35 L= 16)   382    104        1 ( 420/1039- 509/1108) -39.8 E37-E31 N18-N23
     10  745 1044 (N14E16 L= 35)   498    130        0 ( 693/1019- 809/1082) -45.0 E18-E12 N11-N17
     11 1870  981 (S18W59 L=110)  1368    151        2 (1831/ 916-1909/1045) -59.0 W54-W64 S23-S13
     12  668  932 (N09E23 L= 28)   296     68        4 ( 649/ 898- 690/ 972)  76.0 E24-E21 N07-N11
     13 1841  891 (S22W54 L=105)   508     52        1 (1819/ 863-1860/ 926) -51.3 W51-W56 S24-S19
     14 1398  758 (S17W16 L= 67)  1367    435        1 (1248/ 683-1529/ 900)  15.9 W06-W27 S19-S11
     15  348  787 (N06E49 L=  1)   692    153        2 ( 325/ 722- 367/ 873)  77.5 E51-E47 N02-N11
     16 1585  515 (S39W28 L= 78)  4751    636        6 (1338/ 355-1835/ 671)   3.2 W04-W58 S45-S33
     17  685  588 (S11E31 L= 20)   338     83        2 ( 645/ 570- 728/ 611) -26.6 E32-E28 S12-S09
     18  322  452 (S15E70 L=340)   981    114        0 ( 300/ 413- 358/ 503)  63.4 E74-E67 S17-S11
     19 1159  265 (S46E13 L= 38)   315     61        1 (1133/ 244-1192/ 284)   9.5 E15-E09 S47-S44
  

    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