Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  748 1863 (N39E50 L=181)   492    148        0 ( 732/1832- 765/1897) -29.1 E54-E45 N36-N42
      2  574 1539 (N12E41 L=189)  1052    287        0 ( 530/1405- 645/1658) -80.5 E44-E38 N03-N21
      3 1332 1233 (N12W12 L=243)  1980    451        2 (1238/1036-1463/1459) -39.6 W01-W24 N04-N23
      4 1495 1411 (N28W19 L=249)   279     40        0 (1471/1385-1514/1437) -26.6 W16-W20 N26-N29
      5  723 1264 (S01E23 L=208)  2195    654        6 ( 489/1071- 821/1366) -18.4 E40-E13 S10-N04
      6  239 1272 (S09E60 L=171)   367     70        0 ( 216/1245- 262/1304)  63.4 E61-E58 S10-S06
      7 1799  991 (N13W50 L=281)  3007    469        5 (1764/ 792-1830/1205) -60.5 W44-W59 N02-N25
      8  630 1080 (S13E24 L=206)  1007    198        2 ( 541/1033- 731/1130)  45.0 E29-E19 S17-S08
      9  167 1032 (S24E61 L=170)   542    172        1 ( 136/ 974- 207/1087)  54.5 E63-E57 S27-S20
     10 1008  714 (S25W07 L=238)  1136    383        1 ( 825/ 686-1206/ 752)  27.6 E03-W18 S30-S18
     11  471  708 (S39E30 L=201)  1544    354        4 ( 335/ 689- 618/ 735)  14.0 E42-E18 S42-S35
  

    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