Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1585 1518 (N34W45 L= 88)   512     74        2 (1558/1487-1613/1555) -78.7 W43-W47 N31-N36
      2  411 1489 (N34E48 L=355)   915    241        0 ( 353/1442- 476/1532) -21.8 E55-E40 N30-N36
      3  760 1387 (N27E18 L= 25)   232     88        1 ( 713/1370- 805/1400)   0.0 E21-E14 N25-N27
      4 1665 1222 (N14W45 L= 88)   505    155        3 (1627/1146-1700/1294) -66.0 W42-W47 N09-N18
      5  282 1210 (N15E54 L=349)   973    219        3 ( 256/1117- 316/1281) -51.3 E57-E49 N09-N19
      6  951 1180 (N14E04 L= 39)   192     81        0 ( 911/1142- 990/1214) -38.7 E06-E01 N11-N15
      7 1361  973 (N00W21 L= 64)   252     89        1 (1323/ 943-1410/ 996)  18.4 W18-W24 S01-N01
      8  575  839 (S07E29 L= 14)   314    113        1 ( 530/ 804- 630/ 873) -29.7 E32-E25 S09-S05
      9  246  656 (S20E61 L=342)   799    218        1 ( 223/ 598- 267/ 706)  48.8 E66-E58 S24-S16
     10 1457  650 (S21W29 L= 72)   275    120        0 (1400/ 629-1511/ 666)   7.1 W24-W32 S22-S19
     11 1227  571 (S26W13 L= 56)   247     85        0 (1182/ 557-1267/ 589)  -9.5 W09-W15 S26-S24
     12  438  543 (S27E45 L=358)   694    135        1 ( 400/ 514- 484/ 565)  24.0 E49-E40 S29-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