Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1569 1522 (N27W49 L= 84)   533     93        0 (1528/1479-1601/1560) -74.1 W47-W50 N23-N30
      2 1793 1381 (N13W64 L=100)   625    197        0 (1753/1360-1829/1414)  -5.2 W59-W70 N12-N14
      3  532 1308 (N30E30 L=  5)  1719    532        6 ( 373/1259- 691/1367)   2.3 E43-E18 N27-N32
      4  306 1193 (N24E49 L=346)   371     49        1 ( 284/1165- 331/1218)  45.0 E51-E47 N22-N25
      5 1114 1172 (N14W08 L= 43)   505    147        4 (1064/1142-1183/1209) -45.0 W04-W11 N10-N16
      6  509  770 (S03E37 L=359)   814    177        0 ( 465/ 686- 546/ 861) -73.3 E38-E35 S08-N02
      7  755  753 (S07E20 L= 15)   535    167        2 ( 673/ 716- 837/ 779) -21.8 E25-E15 S09-S04
      8 1589  689 (S24W33 L= 69)   819    166        4 (1521/ 640-1648/ 736)  15.3 W27-W38 S25-S21
      9  498  629 (S12E41 L=355)   798    163        3 ( 463/ 565- 540/ 686)  45.0 E45-E37 S15-S08
     10  685  462 (S25E32 L=  3)   376    142        1 ( 653/ 406- 716/ 504)  40.6 E36-E29 S28-S22
     11  558  282 (S38E53 L=343)  1017    174        2 ( 533/ 254- 585/ 316)  68.2 E58-E50 S41-S35
  

    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