Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1409 1738 (N35W57 L=143)  2646    340        2 (1344/1622-1493/1828)  15.7 W47-W80 N27-N39
      2 1745 1062 (S12W47 L=132)   294     64        0 (1729/1041-1769/1094) -63.4 W45-W48 S13-S09
      3 1642  909 (S17W35 L=120)  1028    246        1 (1620/ 796-1662/1024)  61.7 W31-W38 S23-S10
      4  480  971 (N17E35 L= 50)   270     54        0 ( 454/ 943- 505/ 994) -63.4 E36-E33 N14-N18
      5  264  807 (N11E58 L= 27)  2354    397        5 ( 159/ 667- 325/ 967) -59.9 E68-E54 N02-N21
      6 1477  563 (S33W16 L=102)   320     64        0 (1446/ 536-1498/ 592)  36.9 W13-W18 S34-S31
      7 1005  440 (S27E18 L= 67)   284     37        0 ( 986/ 415-1029/ 460)  33.7 E20-E16 S28-S25
      8 1411  339 (S49W06 L= 92)   332     19        0 (1391/ 324-1432/ 355)  14.0 W04-W08 S51-S47
      9 1314  272 (S52E05 L= 80)   388     43        0 (1287/ 254-1342/ 287)   0.0 E08-E02 S53-S50
     10  813  257 (S35E42 L= 43)   479     51        0 ( 782/ 236- 846/ 282)   7.1 E45-E37 S36-S33
  

    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