Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  780 1813 (N56E39 L=308)  1447    109        2 ( 739/1752- 819/1850) -51.3 E46-E32 N50-N60
      2  971 1595 (N39E09 L=338)   247     46        1 ( 954/1572- 997/1611)  33.7 E10-E06 N37-N40
      3  522 1494 (N28E42 L=305)   657    142        1 ( 460/1461- 580/1537) -13.0 E49-E36 N26-N30
      4 1514 1408 (N29W33 L= 20)   682    198        1 (1455/1339-1583/1479) -41.6 W29-W38 N25-N34
      5  331 1390 (N19E56 L=291)  1401    170        2 ( 272/1322- 371/1456) -22.4 E66-E49 N15-N23
      6  787 1323 (N19E18 L=329)   579    227        0 ( 706/1250- 866/1372) -26.6 E23-E11 N15-N21
      7 1360 1074 (N07W21 L=  8)   444     93        2 (1328/1035-1385/1115)  51.3 W18-W22 N04-N09
      8 1275  696 (S17W19 L=  5)   227     69        0 (1254/ 662-1296/ 729)  63.4 W17-W19 S18-S14
      9 1467  354 (S40W46 L= 33)   514     80        0 (1443/ 322-1485/ 395) -60.3 W43-W49 S43-S36
     10 1113  316 (S46W15 L=  2)   331     47        1 (1088/ 298-1137/ 339) -18.4 W11-W17 S47-S44
     11 1266  297 (S47W30 L= 17)   472     26        0 (1242/ 272-1287/ 318) -23.2 W26-W33 S48-S44
  

    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