Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  904 1544 (N26E12 L=196)   913    260        0 ( 799/1474- 996/1592) -13.0 E18-E05 N22-N29
      2 1624 1571 (N36W43 L=251)   641     72        0 (1575/1553-1661/1591)  16.7 W37-W47 N33-N37
      3  259 1391 (N14E59 L=150)  3563    362        7 ( 199/1246- 355/1552) -46.5 E70-E46 N05-N25
      4  359 1425 (N16E49 L=159)   964    170        4 ( 308/1379- 417/1513)  66.0 E53-E46 N13-N22
      5  429 1421 (N16E43 L=165)   488     66        2 ( 398/1388- 465/1451)  -9.5 E46-E40 N14-N18
      6  485 1399 (N15E38 L=170)   522    106        0 ( 443/1359- 520/1444)  59.0 E40-E35 N12-N17
      7 1468 1288 (N14W26 L=234)   335     43        0 (1444/1256-1488/1312) -56.3 W24-W27 N12-N15
      8 1289 1279 (N12W14 L=223)   194     94        0 (1247/1257-1329/1295)  31.0 W11-W16 N10-N13
      9  190 1076 (S05E60 L=148)   801    150        1 ( 132/1035- 248/1117)  18.4 E67-E53 S07-S02
     10  589 1037 (S07E27 L=181)   297     83        1 ( 563/1007- 608/1076) -76.0 E28-E25 S08-S04
     11  184  885 (S16E62 L=146)  3152    305        4 ( 136/ 787- 261/1017)  51.8 E69-E53 S22-S08
     12  507  976 (S10E32 L=176)   267     57        1 ( 486/ 948- 533/1000)  56.3 E33-E30 S12-S08
     13  419  880 (S16E39 L=169)  2680    551        6 ( 245/ 778- 601/ 952)  17.9 E56-E25 S22-S11
     14  200  684 (S29E66 L=142)   701     62        1 ( 183/ 647- 219/ 723)  45.0 E69-E63 S30-S26
     15 1313  518 (S33W25 L=233)   450    104        0 (1251/ 494-1361/ 543)   0.0 W20-W29 S34-S32
     16  983  463 (S40W02 L=210)   252     94        0 ( 935/ 449-1027/ 477)  18.4 E01-W05 S41-S38
     17 1473  244 (S51W61 L=269)   604     53        0 (1452/ 219-1497/ 272) -10.6 W54-W70 S52-S49
     18 1371  243 (S53W48 L=256)   758    108        0 (1321/ 226-1426/ 260)   4.4 W40-W54 S54-S51
  

    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