Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1474 1511 (N28W32 L=227)   256     96        0 (1430/1479-1513/1531) -31.0 W28-W34 N25-N28
      2  277 1274 (N13E53 L=142)  2298    410        6 ( 186/1151- 360/1447)  74.5 E61-E47 N06-N24
      3 1723 1412 (N22W52 L=247)   320     61        1 (1696/1391-1747/1434) -26.6 W49-W54 N20-N23
      4  680 1314 (N15E22 L=173)   674    208        3 ( 592/1292- 785/1339)  -4.4 E27-E14 N13-N16
      5  337 1169 (N07E46 L=149)   778    157        2 ( 290/1098- 383/1228)  53.1 E49-E42 N02-N10
      6  901  839 (S14E08 L=187)   950    311        1 ( 837/ 707- 938/ 980)  70.6 E12-E05 S22-S05
      7 1374  821 (S15W22 L=217)   246     78        0 (1349/ 778-1396/ 847) -53.1 W20-W23 S17-S13
      8 1458  800 (S17W28 L=223)   401    152        1 (1399/ 775-1507/ 827)   0.0 W23-W31 S17-S14
      9  898  472 (S38E10 L=185)   727    207        1 ( 791/ 452- 995/ 504)   7.1 E18-E02 S39-S35
     10 1691  442 (S39W63 L=258)   613     92        0 (1655/ 420-1731/ 468)  -9.5 W55-W73 S40-S37
     11  872  303 (S52E15 L=180)   549    129        0 ( 807/ 283- 930/ 321)  -4.8 E21-E09 S53-S50
  

    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