Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  848 1650 (N51W07 L=290)   522    109        1 ( 809/1612- 884/1691)  24.0 W02-W11 N48-N52
      2 1193 1607 (N37W30 L=313)   515     94        2 (1150/1575-1239/1626)  -7.1 W26-W34 N35-N38
      3 1587 1529 (N20W54 L=336)  1211    177        4 (1507/1481-1661/1555) -28.6 W48-W61 N16-N23
      4 1395 1424 (N20W34 L=317)   399    151        0 (1381/1351-1414/1491)  60.9 W31-W37 N15-N24
      5  320 1287 (N38E48 L=234)  1039    195        2 ( 237/1218- 383/1332) -25.5 E61-E40 N32-N42
      6 1512 1269 (N08W36 L=319)   226     56        0 (1486/1240-1537/1296) -76.0 W34-W36 N06-N10
      7 1867 1181 (S08W65 L=347)   755     66        1 (1841/1152-1889/1215)  31.0 W60-W68 S09-S05
      8  679 1198 (N25E17 L=265)   557    200        0 ( 582/1178- 766/1213) -21.0 E24-E11 N22-N27
      9 1217 1165 (N11W14 L=297)   362    152        0 (1135/1148-1272/1185) -18.4 W09-W18 N08-N12
     10 1525 1057 (S03W30 L=313)   695    145        1 (1488/ 982-1562/1121)  45.0 W26-W34 S07-N00
     11 1054 1071 (N09W03 L=286)   196     63        1 (1026/1051-1082/1096)  11.3 W00-W05 N08-N10
     12  282 1030 (N22E52 L=231)  1112    158        1 ( 227/ 970- 337/1086) -38.7 E58-E47 N18-N26
     13 1847 1026 (S15W57 L=340)   405     68        1 (1824/1000-1880/1047)   0.0 W54-W62 S16-S14
     14 1547  731 (S23W26 L=308)  1385    279        1 (1458/ 610-1620/ 839)  33.7 W17-W32 S29-S18
     15  914  782 (S04E12 L=271)   197     70        0 ( 890/ 755- 940/ 816) -78.7 E12-E10 S06-S01
     16  803  713 (S06E20 L=263)   246     69        0 ( 771/ 697- 838/ 735) -45.0 E21-E18 S07-S03
     17 1340  375 (S43W04 L=287)  2501    351        8 (1217/ 278-1434/ 517)  28.8 E08-W13 S49-S34
     18  808  288 (S34E38 L=245)   542     76        0 ( 771/ 263- 850/ 325)  16.7 E42-E32 S35-S31
     19 1066  203 (S50E26 L=256)   559     89        1 (1027/ 181-1115/ 223)  -7.1 E30-E21 S51-S47
  

    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