Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  385 1483 (N04E55 L=206)  1018    201        0 ( 333/1388- 436/1557)  90.0 E56-E52 S02-N09
      2 1520 1519 (N38W17 L=278)   503     86        1 (1480/1500-1564/1543)  15.9 W14-W21 N36-N39
      3  234 1294 (S10E62 L=199)  1162    197        3 ( 164/1254- 286/1342)  24.8 E70-E56 S12-S06
      4  465 1252 (S08E40 L=221)   706    215        1 ( 434/1162- 532/1317)  74.7 E42-E36 S13-S02
      5 1148 1258 (N10W00 L=261)   215     57        1 (1121/1233-1172/1289) -45.0 E01-W02 N08-N11
      6  312 1030 (S23E47 L=214)   883    208        0 ( 273/ 941- 387/1096)  54.0 E50-E41 S29-S18
      7 1216  767 (S15W18 L=279)   230     57        0 (1192/ 745-1233/ 799) -33.7 W16-W19 S16-S13
      8 1517  660 (S11W40 L=301)   243     62        0 (1494/ 626-1530/ 690) -21.8 W37-W42 S12-S09
      9  452  585 (S48E30 L=231)   686    117        2 ( 395/ 566- 509/ 608)   8.7 E36-E23 S50-S45
     10 1446  547 (S19W41 L=302)   366     87        0 (1400/ 531-1493/ 559)  26.6 W37-W43 S20-S17
     11 1334  491 (S26W37 L=298)   432     91        2 (1285/ 468-1380/ 510)  15.9 W33-W40 S27-S24
     12  992  331 (S47W26 L=286)  1208    225        0 ( 902/ 275-1097/ 383)   2.7 W16-W37 S48-S45
  

    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