Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  830 1648 (N39E25 L=181)   454    117        1 ( 783/1619- 883/1674) -11.3 E29-E19 N37-N40
      2 1582 1484 (N37W39 L=245)  1954    328        5 (1539/1344-1614/1626)  76.6 W34-W45 N27-N48
      3 1045 1558 (N36E06 L=200)   356     50        0 (1022/1528-1066/1585) -53.1 E07-E04 N33-N37
      4  585 1484 (N24E38 L=168)   407    137        0 ( 552/1443- 639/1506)  38.7 E39-E33 N21-N26
      5 1104 1440 (N28W00 L=206)   238     72        1 (1067/1417-1133/1459) -14.0 E02-W02 N26-N29
      6  910 1401 (N23E12 L=194)   257     61        2 ( 877/1379- 945/1421)   0.0 E14-E09 N21-N24
      7 1244 1211 (N15W12 L=217)   299    137        0 (1211/1167-1298/1262) -39.8 W09-W15 N12-N17
      8  359  930 (S13E45 L=161)   592    104        1 ( 326/ 894- 389/ 961) -45.0 E47-E41 S14-S10
      9 1505  895 (S02W32 L=238)   339    129        0 (1450/ 854-1544/ 940) -20.6 W27-W35 S03-N00
     10 1348  754 (S12W24 L=230)   936    303        0 (1222/ 684-1450/ 839)  42.7 W16-W29 S17-S05
     11 1318  774 (S11W22 L=227)   214    111        0 (1269/ 749-1342/ 814)  45.0 W18-W22 S13-S08
     12 1128  795 (S12W09 L=215)   194     44        0 (1101/ 779-1154/ 811)  33.7 W07-W10 S12-S10
     13  527  698 (S26E31 L=175)   250     48        0 ( 509/ 679- 548/ 722) -45.0 E32-E28 S26-S23
     14 1184  589 (S24W16 L=222)   475     99        1 (1130/ 554-1218/ 617)  -8.1 W12-W19 S26-S22
     15  452  526 (S39E41 L=165)   325     37        1 ( 437/ 505- 472/ 548)  45.0 E42-E38 S40-S37
     16  718  535 (S34E16 L=190)   301     80        1 ( 675/ 522- 759/ 548)  15.9 E19-E12 S35-S33
     17  378  494 (S43E52 L=154)  1355    126        3 ( 310/ 472- 427/ 536)   8.5 E65-E45 S45-S39
     18 1095  371 (S42W16 L=222)   304     69        1 (1048/ 360-1128/ 386)   9.5 W12-W18 S43-S40
     19 1045  233 (S56W18 L=224)   391     32        1 (1022/ 211-1067/ 254)  45.0 W15-W20 S57-S52
  

    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