Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  821 1870 (N68W22 L= 27)  1327    117        2 ( 787/1820- 864/1921)  17.2 W13-W42 N62-N72
      2  909 1876 (N64W32 L= 37)   422     53        0 ( 895/1846- 925/1907)  30.3 W27-W40 N60-N67
      3 1162 1752 (N44W36 L= 41)   763    132        1 (1129/1705-1194/1806)  28.6 W30-W41 N40-N47
      4  923 1745 (N51W18 L= 23)   467     69        1 ( 882/1725- 964/1760) -15.9 W14-W21 N49-N53
      5 1448 1639 (N28W48 L= 53)   676    124        1 (1397/1612-1509/1676) -54.5 W45-W52 N24-N31
      6 1282 1563 (N28W32 L= 37)   779    324        2 (1162/1518-1416/1644) -47.3 W26-W39 N22-N35
      7 1564 1430 (N13W45 L= 50)   627    141        1 (1527/1362-1599/1494)  63.4 W43-W47 N08-N16
      8 1618 1430 (N11W49 L= 54)   277     33        0 (1601/1407-1640/1450)  11.3 W46-W51 N10-N12
      9  915 1417 (N29W03 L=  9)   580    209        0 ( 827/1348- 971/1447) -60.9 E01-W06 N23-N32
     10 1187 1421 (N22W20 L= 25)   474    239        1 (1074/1406-1295/1435) -19.7 W12-W26 N18-N25
     11 1122 1360 (N20W14 L= 19)   270     96        0 (1075/1346-1171/1378) -38.7 W11-W16 N18-N22
     12  343 1122 (N24E43 L=322)   475     85        0 ( 310/1085- 384/1158) -45.0 E45-E39 N21-N27
     13  308 1001 (N18E48 L=317)   614     58        2 ( 277/ 968- 342/1033) -38.7 E50-E44 N15-N20
     14  204  906 (N14E60 L=305)   718    162        1 ( 189/ 833- 221/ 988)  66.0 E64-E58 N10-N19
     15  831  670 (S12E20 L=345)   293    114        1 ( 805/ 606- 859/ 701)  39.8 E23-E17 S15-S09
     16 1708  649 (S37W40 L= 45)   509     55        1 (1683/ 623-1734/ 678)  26.6 W37-W42 S38-S35
     17 1285  442 (S38W01 L=  7)  1605    348        5 (1163/ 326-1392/ 552)  22.2 E11-W11 S44-S32
     18 1455  270 (S61W14 L= 19)   719     81        1 (1415/ 249-1498/ 293)  -4.8 W08-W20 S63-S58
     19 1321  186 (S65E07 L=358)   515     69        0 (1290/ 169-1349/ 207)  -7.1 E10-E01 S67-S62
  

    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