Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1241 1627 (N43W06 L=318)   615     99        2 (1215/1580-1269/1670) -66.8 W03-W08 N39-N46
      2 1783 1299 (N29W55 L=  8)  1526    314        1 (1698/1186-1849/1411) -40.9 W47-W62 N22-N35
      3 1007 1341 (N19E06 L=307)   202     97        0 ( 968/1308-1047/1367)  51.3 E07-E03 N16-N21
      4 1754 1256 (N25W51 L=  4)   534    143        0 (1710/1198-1786/1308) -40.6 W46-W54 N22-N28
      5 1214 1220 (N15W09 L=322)   381     99        0 (1175/1184-1252/1261) -31.0 W05-W11 N12-N16
      6 1343  824 (S07W23 L=336)   227    120        0 (1286/ 803-1387/ 850)  -7.1 W18-W26 S07-S06
      7 1360  686 (S15W27 L=340)   869    246        0 (1290/ 611-1397/ 745) -83.7 W21-W29 S20-S11
      8 1467  674 (S15W35 L=348)   522    119        3 (1433/ 619-1499/ 724)  90.0 W32-W36 S18-S11
      9 1512  607 (S18W40 L=353)   838    258        2 (1408/ 547-1602/ 670)  36.9 W33-W45 S23-S12
     10  735  608 (S31E13 L=300)  1388    427        1 ( 569/ 554- 925/ 666)  26.6 E26-W00 S37-S23
     11 1597  393 (S31W62 L= 15)  3717    280        5 (1476/ 347-1696/ 465)  26.6 W49-W74 S34-S24
     12 1081  344 (S44W19 L=332)  1339    226        4 ( 997/ 319-1200/ 372)  26.6 W12-W27 S47-S40
  

    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