Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1703 1366 (N11W54 L= 54)   763    198        2 (1650/1291-1743/1438) -66.0 W50-W57 N06-N15
      2 1535 1307 (N12W38 L= 38)   256     60        1 (1507/1279-1558/1327) -45.0 W36-W40 N10-N13
      3 1331 1274 (N15W24 L= 24)   218     51        0 (1304/1250-1354/1296) -56.3 W22-W25 N13-N16
      4 1349 1230 (N12W24 L= 23)   256    123        0 (1298/1199-1414/1252) -39.8 W21-W27 N09-N14
      5  184 1157 (N27E61 L=299)   677    132        2 ( 162/1124- 214/1197) -63.4 E63-E57 N25-N29
      6  771 1109 (N16E14 L=346)   365    133        1 ( 706/1071- 836/1142)   0.0 E18-E09 N14-N17
      7  544 1101 (N19E29 L=330)   568    164        3 ( 477/1058- 606/1131) -36.9 E33-E24 N16-N22
      8  647  988 (N11E24 L=336)   260     71        3 ( 627/ 948- 666/1025) -90.0 E24-E22 N08-N13
      9  665  825 (N01E25 L=334)   860    290        3 ( 546/ 716- 793/ 899)  11.9 E35-E16 S03-N04
     10  686  622 (S11E29 L=331)   325    157        1 ( 619/ 602- 751/ 654)   0.0 E33-E23 S12-S09
     11  334  480 (S17E64 L=295)  1215    283        2 ( 308/ 395- 355/ 567)  65.2 E69-E60 S23-S10
  

    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