Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1541 1539 (N29W49 L= 62)   723    224        2 (1478/1463-1612/1581)  21.8 W39-W54 N25-N31
      2 1171 1446 (N30W17 L= 30)   477    239        1 (1151/1349-1195/1547)  72.9 W14-W19 N23-N36
      3 1661 1461 (N22W55 L= 69)   364     61        1 (1636/1433-1688/1484) -45.0 W53-W57 N19-N23
      4  442 1287 (N29E39 L=334)  1546    357        3 ( 406/1141- 512/1441) -64.5 E44-E32 N18-N39
      5 1654 1388 (N18W51 L= 64)   317     87        1 (1619/1357-1679/1422) -59.0 W49-W52 N15-N20
      6 1548  990 (S04W33 L= 46)   368    142        0 (1525/ 923-1583/1039)  63.4 W31-W36 S07-S01
      7 1754  862 (S15W49 L= 62)   566     69        1 (1728/ 837-1780/ 891)  21.8 W46-W51 S17-S14
      8  654  557 (S19E31 L=342)   255     72        0 ( 625/ 536- 692/ 578) -26.6 E33-E28 S20-S17
      9  586  510 (S21E38 L=335)   294     83        1 ( 561/ 481- 625/ 524)   0.0 E40-E34 S23-S20
     10 1010  495 (S27E09 L=  5)   246    109        0 ( 961/ 485-1059/ 516)   7.1 E12-E04 S28-S26
     11  614  386 (S31E42 L=331)   326     63        0 ( 594/ 357- 637/ 417) -90.0 E43-E40 S33-S28
  

    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