Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  452 1605 (N36E56 L=130)   643    114        1 ( 424/1558- 486/1652) -24.8 E63-E49 N32-N39
      2 1362 1460 (N32W22 L=208)  1610    470        4 (1233/1317-1538/1561) -33.7 W12-W34 N23-N39
      3  662 1406 (N24E28 L=158)   572    186        0 ( 624/1333- 712/1485) -42.0 E32-E22 N19-N29
      4 1669 1259 (N20W45 L=231)   445    103        2 (1623/1239-1705/1297)  -8.1 W41-W48 N18-N22
      5 1725  826 (S07W51 L=236)  1092    181        2 (1648/ 760-1790/ 872)  32.5 W44-W57 S11-S03
      6 1312  723 (S15W20 L=206)   304    104        1 (1293/ 683-1331/ 771)  71.6 W19-W21 S17-S11
      7 1472  634 (S20W33 L=219)   398    111        0 (1432/ 601-1503/ 678)  59.0 W30-W35 S22-S16
      8  479  550 (S32E39 L=146)   707    152        2 ( 415/ 536- 547/ 572)   8.7 E45-E32 S33-S30
  

    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