Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1071 1337 (N12E03 L=144)   705    246        2 (1042/1224-1098/1457) -66.8 E06-W00 N05-N19
      2  357 1423 (N06E53 L= 95)   352     85        0 ( 317/1401- 391/1453)  31.0 E55-E50 N04-N08
      3  347 1359 (N02E51 L= 96)   515    111        2 ( 307/1326- 380/1382)  -8.1 E55-E48 N00-N03
      4 1712 1221 (N20W40 L=188)   720    198        4 (1687/1131-1743/1287)  79.7 W38-W42 N13-N24
      5 1337 1167 (N07W15 L=163)   300    121        0 (1308/1115-1367/1225)  81.9 W14-W16 N04-N11
      6 1715  923 (N03W46 L=193)   550    193        1 (1672/ 874-1768/ 998) -20.0 W40-W51 N00-N05
      7 1748  714 (S08W56 L=203)   521    110        2 (1719/ 661-1783/ 754)  66.8 W52-W57 S11-S04
      8 1195  674 (S23W18 L=166)   971    338        2 (1044/ 645-1337/ 746)  36.3 W09-W24 S27-S16
      9 1581  550 (S21W48 L=196)   640    164        3 (1505/ 516-1650/ 583)  41.2 W43-W52 S24-S17
     10 1398  545 (S26W35 L=183)   282     64        0 (1365/ 523-1429/ 562)  53.1 W33-W37 S28-S24
     11 1465  276 (S40W56 L=204)  1923    366        5 (1356/ 254-1562/ 302)  19.4 W46-W64 S43-S36
  

    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