Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  748 1764 (N30E44 L= 60)   411     68        1 ( 711/1741- 779/1785)   0.0 E47-E40 N28-N31
      2 1023 1719 (N35E23 L= 81)   439     93        0 ( 988/1677-1049/1763) -32.0 E27-E19 N32-N37
      3 1088 1400 (N16E07 L= 97)   551    163        0 (1049/1329-1119/1480) -41.2 E10-E02 N12-N19
      4 1214 1358 (N17W02 L=106)   473    133        0 (1168/1334-1280/1396)  18.4 E01-W05 N15-N18
      5 1320 1260 (N14W10 L=114)   288    102        0 (1302/1214-1348/1306) -38.7 W08-W13 N11-N16
      6 1737  777 (N01W52 L=156)   766     86        1 (1711/ 737-1761/ 822)  -8.1 W48-W55 S01-N02
      7 1711  650 (S07W56 L=160)   384     36        0 (1693/ 626-1730/ 676)  90.0 W54-W57 S08-S05
      8 1045  535 (S34W17 L=121)  2078    486        2 ( 865/ 480-1264/ 651)  49.4 W06-W25 S42-S20
      9 1449  508 (S23W44 L=148)   507     73        2 (1419/ 479-1479/ 537)   0.0 W40-W47 S24-S20
     10 1475  432 (S26W51 L=155)   309     71        0 (1443/ 399-1499/ 460)  90.0 W49-W52 S28-S23
     11 1175  392 (S38W33 L=137)   289     70        0 (1161/ 351-1191/ 426) -51.3 W31-W35 S40-S35
     12 1108  303 (S46W36 L=140)   405     98        0 (1062/ 292-1164/ 317)  29.7 W32-W39 S48-S43
  

    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