Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1816 (N43W59 L= 78)   629     55        1 (1278/1780-1326/1840)  26.6 W54-W63 N40-N45
      2 1292 1465 (N24W30 L= 49)   406     69        1 (1258/1441-1316/1497)   0.0 W26-W32 N21-N25
      3 1820 1330 (S01W64 L= 84)   985     94        1 (1772/1309-1879/1349) -10.0 W58-W75 S02-N01
      4 1038 1222 (N17W06 L= 25)   874    263        1 ( 983/1117-1125/1326)  26.6 W00-W14 N12-N20
      5  832 1246 (N24E05 L= 14)   287    112        0 ( 784/1225- 890/1273) -45.0 E07-E02 N21-N26
      6 1414 1241 (N08W28 L= 48)   215     49        0 (1398/1219-1431/1265)  26.6 W27-W30 N06-N09
      7  674 1142 (N22E18 L=  1)   605    213        0 ( 574/1089- 747/1175) -38.7 E24-E14 N17-N25
      8  415 1089 (N25E38 L=341)   373    111        1 ( 395/1035- 443/1136) -74.1 E39-E35 N21-N28
      9  176 1038 (N26E65 L=315)   575     75        1 ( 153/ 987- 194/1075) -59.0 E68-E62 N23-N28
     10 1739  951 (S18W43 L= 63)   858    204        3 (1650/ 927-1843/ 978) -24.0 W36-W54 S22-S13
     11 1372  898 (S10W16 L= 35)   367     97        2 (1343/ 849-1403/ 945)  38.7 W13-W19 S11-S07
     12 1683  901 (S19W38 L= 57)   315     85        1 (1651/ 872-1726/ 930)  -8.1 W34-W41 S20-S18
     13  231  780 (N10E62 L=317)   616     86        0 ( 210/ 742- 260/ 830) -80.5 E64-E59 N07-N13
     14 1252  655 (S20W03 L= 22)   378     90        2 (1233/ 615-1272/ 704) -90.0 W01-W04 S23-S17
     15  893  584 (S15E20 L=360)   588    191        0 ( 857/ 495- 937/ 656)  38.7 E24-E14 S19-S11
     16  660  428 (S19E41 L=338)   491    124        0 ( 639/ 379- 687/ 474)  56.3 E43-E37 S22-S16
     17  955  383 (S30E24 L=355)   739    127        1 ( 910/ 334-1004/ 423)  16.7 E29-E19 S32-S27
     18  599  359 (S22E50 L=329)   815     87        2 ( 561/ 328- 627/ 396)  18.4 E54-E45 S24-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