Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  905 1716 (N41E27 L=191)   512     71        1 ( 882/1678- 927/1755) -35.5 E31-E23 N37-N43
      2 1293 1432 (N29W09 L=227)  2037    598        1 (1110/1186-1443/1577) -31.8 E06-W23 N16-N35
      3 1195 1191 (N12W07 L=225)   625    175        1 (1128/1161-1284/1225)   0.0 W02-W12 N09-N13
      4 1787 1050 (N15W53 L=271)  1097    193        2 (1759/ 971-1832/1117) -48.8 W49-W59 N11-N19
      5  267  938 (S20E53 L=165)  2906    411        2 ( 143/ 795- 381/1077)  33.3 E75-E42 S31-S09
      6 1792  869 (N04W57 L=275)   487     75        0 (1768/ 833-1815/ 906)  59.0 W55-W58 N02-N07
      7 1121  717 (S18W12 L=230)   342     44        3 (1094/ 692-1150/ 741) -33.7 W10-W13 S19-S16
      8 1610  646 (S12W47 L=265)  1554    360        4 (1457/ 575-1729/ 696)  30.1 W38-W59 S19-S08
      9 1586  511 (S21W52 L=270)   465    117        1 (1559/ 475-1633/ 541)  45.0 W49-W55 S23-S18
     10  857  293 (S54W06 L=224)   688     91        2 ( 803/ 265- 903/ 325)  -5.2 W00-W11 S56-S52
     11  706  282 (S61E10 L=208)   405     45        0 ( 677/ 269- 736/ 297)  20.6 E14-E06 S62-S58
  

    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