Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1325 1546 (N33W03 L= 29)   232     51        1 (1295/1529-1353/1566)  14.0 W00-W04 N31-N34
      2 1522 1237 (N20W23 L= 49)   716    182        0 (1458/1157-1560/1317) -31.0 W17-W27 N16-N22
      3 1853 1022 (N20W55 L= 81)  1494    219        5 (1827/ 928-1878/1144) -59.7 W50-W59 N14-N26
      4  286 1085 (S21E50 L=336)   352     80        0 ( 275/1045- 298/1132) -78.7 E51-E48 S23-S18
      5 1192 1005 (S03W10 L= 36)   286     81        1 (1172/ 959-1208/1052) -68.2 W07-W10 S05-N00
      6  683  755 (S31E13 L= 13)  1452    439        1 ( 486/ 700- 887/ 830)  10.1 E28-W00 S33-S25
      7 1424  580 (S18W38 L= 64)  1199    339        2 (1323/ 487-1550/ 731)  81.5 W34-W41 S26-S06
      8 1145  475 (S33W26 L= 52)   232     81        0 (1109/ 460-1192/ 497)   9.5 W23-W29 S33-S32
  

    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