Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1479 1692 (N30W54 L=  6)   797    209        0 (1385/1645-1570/1722) -42.0 W47-W58 N25-N34
      2  412 1632 (N58E44 L=268)   873    140        1 ( 352/1587- 462/1675)  18.4 E53-E36 N55-N61
      3 1459 1642 (N28W48 L=  0)   345    109        1 (1410/1622-1503/1659) -14.0 W44-W52 N26-N29
      4 1111 1352 (N20W13 L=325)   380    105        1 (1053/1338-1161/1372) -26.6 W09-W15 N18-N21
      5 1214 1256 (N12W16 L=328)   459     95        1 (1161/1228-1264/1286)  -8.1 W12-W19 N11-N12
      6 1657 1059 (S10W40 L=351)   248     43        0 (1638/1036-1674/1081)  45.0 W38-W41 S10-S08
      7 1752  605 (S41W47 L=359)   543     76        0 (1716/ 587-1784/ 631)   8.1 W41-W51 S43-S39
      8  410  523 (S13E54 L=258)   325     64        0 ( 393/ 500- 432/ 557)  76.0 E55-E52 S14-S10
      9  752  329 (S32E39 L=273)  1411    422        2 ( 623/ 234- 893/ 461)  15.9 E57-E22 S37-S26
     10 1023  311 (S41E22 L=290)   356    101        1 (1006/ 262-1044/ 361)  41.2 E25-E17 S44-S37
  

    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