Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  601 1570 (N24E42 L=168)  3171    544        4 ( 453/1453- 825/1730) -12.4 E71-E21 N19-N32
      2  972 1436 (N22E12 L=198)   278     91        0 ( 935/1420-1021/1456)   0.0 E14-E08 N21-N23
      3  820 1267 (N09E17 L=193)   747    347        4 ( 715/1218- 939/1351)  -7.6 E24-E09 N06-N12
      4  660 1192 (N02E25 L=185)   475    105        2 ( 606/1167- 707/1212)  18.4 E28-E22 N00-N03
      5 1112 1114 (N05W03 L=213)   410    208        3 (1087/1019-1136/1212) -79.7 W02-W04 N00-N11
      6  540  827 (S22E28 L=182)  1005    303        2 ( 492/ 730- 630/ 938) -47.7 E31-E20 S26-S15
      7 1623  455 (S23W60 L=270)   836    134        1 (1557/ 428-1688/ 489)  31.0 W54-W65 S26-S19
      8 1284  369 (S37W37 L=247)   358     43        1 (1261/ 341-1311/ 398) -23.2 W33-W40 S38-S34
      9 1427  373 (S33W48 L=258)   579     82        3 (1389/ 352-1468/ 397)  23.2 W44-W51 S34-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