Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  369 1610 (N29E56 L=322)   810    151        1 ( 352/1557- 383/1664) -69.4 E58-E53 N24-N32
      2 1520 1534 (N32W31 L= 48)   512     91        1 (1494/1507-1547/1570)  78.7 W28-W33 N29-N34
      3  360 1491 (N20E51 L=326)   404     89        0 ( 328/1474- 396/1506)   9.5 E54-E48 N18-N21
      4 1263 1355 (N17W12 L= 30)   452    150        2 (1231/1285-1288/1432) -77.5 W09-W13 N12-N21
      5  631 1358 (N12E27 L=350)   662    142        1 ( 575/1311- 679/1397) -32.0 E31-E23 N09-N14
      6  250 1244 (N04E56 L=321)  1274    187        3 ( 209/1171- 300/1295)  40.6 E59-E51 N00-N07
      7 1514 1185 (N09W29 L= 46)   531    108        0 (1464/1120-1562/1232) -40.6 W25-W32 N05-N11
      8 1919 1069 (N08W67 L= 84)   610     90        0 (1900/1030-1940/1104) -45.0 W64-W70 N05-N09
      9 1791  978 (N00W53 L= 70)   615    128        0 (1771/ 909-1807/1044) -50.2 W50-W55 S03-N03
     10  305  749 (S26E50 L=328)  1887    488        2 ( 232/ 595- 409/ 916)  46.5 E60-E39 S35-S15
     11  483  443 (S46E41 L=336)   330     82        0 ( 465/ 415- 503/ 474)  38.7 E43-E38 S48-S43
  

    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