Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1386 1447 (N31W25 L=354)   555    205        0 (1299/1390-1462/1484) -24.4 W18-W30 N26-N33
      2 1681 1405 (N28W51 L= 20)   509    103        2 (1647/1367-1720/1442) -51.3 W47-W54 N25-N30
      3 1684 1180 (N13W46 L= 15)   484    112        1 (1637/1156-1740/1209)  11.3 W41-W51 N11-N15
      4 1581  846 (S07W37 L=  6)   282     97        0 (1563/ 806-1604/ 894) -59.0 W35-W39 S09-S04
      5 1503  705 (S16W33 L=  2)   416     98        4 (1469/ 674-1538/ 745) -59.0 W30-W35 S18-S13
      6 1714  632 (S21W54 L= 23)   438    150        1 (1667/ 575-1750/ 694)  63.4 W51-W56 S25-S17
      7 1174  676 (S19W10 L=339)   255     68        0 (1142/ 655-1214/ 693)  21.8 W08-W13 S19-S17
      8 1445  622 (S22W30 L=359)   789    164        4 (1385/ 563-1515/ 657)  -5.7 W25-W35 S25-S19
      9  324  613 (S26E55 L=274)   644     96        2 ( 296/ 569- 357/ 652)  28.6 E60-E49 S29-S23
     10  515  537 (S30E38 L=291)  1753    315        1 ( 386/ 493- 654/ 606)  19.7 E52-E24 S34-S24
     11 1273  514 (S30W19 L=348)   594    212        1 (1177/ 462-1343/ 563)  36.0 W12-W23 S33-S25
  

    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