Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1462 1691 (N35W52 L=168)   762    155        1 (1391/1667-1513/1712) -16.7 W46-W56 N32-N37
      2  886 1532 (N41W03 L=119)   648     87        2 ( 849/1499- 928/1566) -51.3 W00-W06 N38-N43
      3 1071 1382 (N27W11 L=127)   243     42        1 (1047/1358-1094/1403) -71.6 W10-W13 N25-N28
      4  638 1285 (N30E20 L= 96)  1053    258        1 ( 544/1208- 727/1376) -51.3 E26-E14 N23-N38
      5 1240 1308 (N19W20 L=136)   752    239        0 (1128/1275-1363/1350) -31.6 W13-W26 N14-N23
      6  355 1258 (N34E43 L= 73)   297     88        0 ( 324/1228- 380/1278) -38.7 E45-E40 N31-N35
      7 1492 1125 (N02W30 L=146)  1364    325        5 (1367/1001-1564/1253) -70.6 W25-W33 S05-N12
      8 1919 1095 (S12W69 L=185)   667     45        0 (1904/1070-1932/1125)  26.6 W66-W71 S14-S10
      9  397 1017 (N18E41 L= 75)   341    105        0 ( 371/ 973- 429/1073) -71.6 E42-E39 N15-N21
     10 1864  998 (S16W58 L=174)   454    104        1 (1845/ 962-1886/1040) -80.5 W57-W61 S18-S12
     11 1507  795 (S17W25 L=141)   663    133        1 (1457/ 737-1543/ 838) -74.1 W21-W26 S21-S14
     12  707  574 (S13E30 L= 86)   511     95        1 ( 682/ 524- 734/ 618)  39.8 E32-E26 S16-S10
     13  416  494 (S14E55 L= 61)  1764    325        2 ( 388/ 394- 443/ 605)  49.4 E61-E48 S21-S07
     14 1202  396 (S37E04 L=112)  1331    269        1 (1135/ 330-1242/ 478)  36.0 E10-W01 S41-S31
     15  749  364 (S28E35 L= 81)   708    156        0 ( 672/ 339- 801/ 398) -20.6 E39-E30 S30-S25
     16  548  364 (S25E50 L= 66)   388     84        0 ( 524/ 347- 584/ 380)   0.0 E53-E47 S26-S23
  

    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