Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1855 (N62W29 L=114)   395     48        0 ( 870/1837- 922/1871)  -8.1 W25-W33 N60-N64
      2 1207 1841 (N48W51 L=135)   444     60        0 (1186/1814-1238/1860)   5.7 W45-W56 N46-N49
      3 1095 1820 (N51W40 L=124)   402     55        0 (1072/1795-1125/1843)   6.3 W36-W45 N49-N52
      4 1512 1750 (N31W65 L=149)  1524    203        1 (1437/1669-1577/1819)  67.4 W60-W70 N25-N37
      5  703 1436 (N36E10 L= 75)   331     95        1 ( 658/1405- 746/1460) -45.0 E12-E07 N33-N38
      6  359 1424 (N45E42 L= 42)  1073    156        2 ( 290/1393- 424/1450) -13.2 E51-E34 N41-N47
      7  960 1430 (N29W07 L= 91)   202    109        1 ( 906/1403- 998/1445)   0.0 W03-W09 N27-N29
      8  623 1295 (N30E18 L= 66)   217     85        0 ( 605/1264- 657/1319)  26.6 E19-E15 N27-N30
      9  293 1180 (N30E48 L= 37)   340     42        1 ( 273/1156- 314/1203)  26.6 E49-E45 N28-N31
     10  333 1094 (N24E44 L= 40)   639    160        2 ( 275/1071- 402/1150) -35.0 E49-E38 N21-N28
     11  424  960 (N14E39 L= 46)   838    198        1 ( 389/ 880- 483/1032)  41.2 E42-E32 N10-N17
     12 1334  946 (S08W15 L=100)   414    117        1 (1312/ 889-1365/ 997)  45.0 W13-W18 S10-S05
     13 1303  682 (S22W08 L= 92)   434    133        1 (1290/ 613-1319/ 747)  76.0 W06-W09 S26-S18
     14 1858  700 (S40W60 L=144)   708     65        0 (1824/ 687-1894/ 714) -14.9 W53-W68 S41-S37
     15 1550  613 (S34W24 L=108)   260     38        0 (1529/ 589-1567/ 636) -63.4 W22-W25 S35-S31
     16  522  307 (S27E61 L= 23)  2214    239        3 ( 463/ 253- 588/ 359)   9.8 E80-E50 S30-S23
     17 1320  286 (S52E03 L= 82)   625     81        1 (1274/ 267-1358/ 305) -33.7 E06-W01 S54-S49
     18 1218  195 (S57E18 L= 66)   627     63        1 (1186/ 172-1247/ 218) -18.4 E20-E13 S60-S54
  

    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